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Abstract. In this paper, we study the controUabihty of the two-dimensional relativistic 
Vlasov-Maxwell system in a torus, by means of an interior control. We give two types 
of results. With the geometric control condition on the control set, we prove the local 
exact controllability of the system in large time. Our proof in this case is based on the 
return method, on some results on the control of the Maxwell equations, and on a suitable 
approximation scheme to solve the non-linear Vlasov-Maxwell system on the torus with an 
absorption procedure. 

Without geometric control condition, but assuming that a strip of the torus is contained 
in the control set and under certain additional conditions on the initial data, we establish 
a controllability result on the distribution function only, also in large time. Here, we need 
some additional arguments based on the asymptotics of the Vlasov-Maxwell system with 
large speed of light and on our previous results concerning the controllability of the Vlasov- 
Poisson system with an external magnetic field [14) . 



1. Introduction 

1.1. The relativistic Vlasov-Maxwell system. Consider the relativistic Vlasov-Maxweh 
system: 



(1.1) 



1 



dtf + v- V^./ + \E+--vAB] ■ V^f = l^G, t>0, X := MVZ^ v G R-\ 



dtE + ccnrlB 



fvdv, dtB + c curl E = 0, 



divE 



fdv 



, f]t=o — fo, E\t=o 
where c > is the speed of light and 



/ dv dx, 

R3xT3 

Eq, Bu 



>|t=o 



div B = 0, 



Bn 



1 + 



is the relativistic velocity. The so-called distribution function f{t, x, v) describes the statistical 
distribution of a population of electrons in a collisionless plasma: the quantity f{t, x, v) dx dv 
can be interpreted as the number of particles at time t whose position is close to x and velocity 
close to V. The term lujG is a source in the Vlasov equation. As usual, E and B stand for the 
electric and magnetic fields, which are solutions of the Maxwell equations. In the equation on 
divE, the term /j^Sxts f dv dx stands for a neutralizing background of fixed ions. We take as 
a convention here that the Lebesgue measure of the torus is equal to 1. 
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In this paper, we wih focus on the two-dimensional version of the Vlasov-Maxwell system, 
which means that {x,v) G T2 X R"^. This corresponds to the physical situation in the usual 
three-dimensional geometry where: 

E = iEiit,x),E2it,x),0), 

^ ' ' B = {0,0, b{t,x)). 

In this framework, noticing that the divergence-free condition is automatically satisfied for B, 
the Vlasov-Maxwell system can be rewritten in the simplified form: 

dtf + v- V^/ +(e + -vH] ■ V,,/ = l^G, t>0,xeT'^,v€ M."^, 



(1.3) 



dtEi - cd^^b = - fvi dv, dtE2 + cd^J) = - f 
dtb + c[d^,E2-d^,Ei] = 0, 
da^^Ei + dx2E2 = / fdv- fdvdx. 

yR2xT2 



V2 dv, 



If V = {vi,V2), we set v^ := ('O2,— ^i)- For simplicity, we shall also denote in the sequel 
cmlb:=VH={d^,b,-d,,b). 

We are interested in the controllability properties of this system, by means of an interior 
control G, localized in the space variable in an open subset co of the torus, but without 
localization in the velocity variable. The basic controllability question is the following one. 
Given (/q, Eq, Bq), (/i, Ei, Bi) belonging to a suitable function space and given some control 
time T > 0, is it possible to steer the system from the first state to the second one in time T 
with a well chosen control function G ? In other words, we ask the solution of (|1.3p to satisfy: 

(1-4) f\t=T = fi, E\t=T = El, B\t=x = Bi. 

Note that in order to preserve electroneutrality, the control G has to satisfy 

luj{x)G{t,x,v)dxdv = 0, yt£[0,T]. 

T2xR2 

1.2. Geometric assumptions and main results of the paper. We prove two kinds of 
controllability results: for the first one we assume that the celebrated geometric control con- 
dition is satisfied for co. For the second one, we assume another type of condition on to, as 
well as other conditions on the initial magnetic field and/or the speed of light. 

1.2.1. With the geometric control condition. The Vlasov-Maxwell can be seen as a hyperbolic 
system, made of a transport equation (the Vlasov part) and a coupling between two wave 
equations (the Maxwell part). Concerning the problem of controllability of the linear wave 
equation, it is well-known that the geometric properties of the control set u play a key role. 
In the seminal works |25| and OH], J. Ranch and M. Taylor and C. Bardos, G. Lebeau and J. 
Ranch introduced the so-called geometric control condition on cj, which roughly means that 
every ray of geometric optics meets co before some fixed time, and proved estimates which 
imply the exact controllability of the wave equation. 

In our case, the geometric control condition can be formulated as follows: 

Definition 1.1. The open subset oj o/T^ satisfies the geometric control condition (GCC) if: 
(1.5) For any x G and any direction e G S^, there exists y G such that x + ye ^ uj. 
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Conversely, it is well known (see N. Burq and P. Gerard [6]) that the geometric control con- 
dition is not only sufficient but is also actually a necessary condition to get exact controllability 
of the wave equation. 

It is hence natural to start by assuming that cj does satisfy the geometric control condi- 
tion. Our first result asserts that with this assumption, local controllability near of the 
Vlasov-Maxwell system holds when the controllability time is large enough. This restriction 
is expected because of the finite propagation speed both for the electromagnetic waves and 
for the charged particles. 

Theorem 1 (With geometric control condition). Assume thatuj satisfies the geometric control 
condition. There exists Tq > 0, such that for any T > Tq, there exists k > satisfying the 
following. If for i = 0,1, {fi, Ei, Bi) € H^{T^ x R'^) x H'-^{T'^) x H^{T'^) are such that fi is 
compactly supported in v and satisfy the compatibility conditions 

(1.6) divE'j = / fidv- fidvdx, 

(1.7) / fodvdx= fidvdx, 
as well as the smallness assumption 

(1.8) \\{fi,Ei,Bi)\\H^<K, f = 0,l, 

then there exists a control G G i7^([0, T] xT^ xM^), supported in to, such that the corresponding 
solution (f,E,B) to the Vlasov-Maxwell system ()1.3p with initial data (fo, Eq, Bq) satisfies: 

f\t = T = /ll -£'|i = T = B^f^rp = Bi. 

Remark 1.1. The solution that we build involves an electric field E and a magnetic field B 
which belong to the Sobolev space C{[0,T]; H^{T'^)), which is imbedded in C{[0,T];C^{T'^)), 
so that characteristics are well defined in the classical sense. 

Remark 1.2. It is possible to generalize this theorem in two directions: relax the smallness 
assumption on Bi and replace the assumption of compact support in v by an assumption of 
exponential decrease in v. We refer to Subsections 15.71 and 15.81 

1.2.2. With the strip assumption. Since GCC is necessary and sufficient for the controllability 
of the wave equation, it seems very involved (although maybe not impossible) to get the result 
of Theorem [1] without assuming GCC. 

Nevertheless, the situation is more favorable if we focus on the system of particles only and 
"forget" about the electromagnetic field: this corresponds to the control of the distribution 
function / only. Indeed, for the kinetic transport equation, we do have more fiexibility than 
for the wave equation and it is natural to try to relax as much as possible the assumptions 
on u. In this work, we are able to prove some results in this direction. Namely, we prove a 
controllability result on the distribution function only with a geometric condition which does 
not involve GCC (and which is not implied by GCC). 

This geometric assumption on the control set will be referred to as the strip assumption 
(loosely speaking, as the name suggests, this corresponds to the case where uj contains a strip). 

Definition 1.2. We say that u satisfies the strip assumption if there exists a straight line of 
whose image Ti by the canonical surjection s : — )■ is closed and included in lo. We 



4 



OLIVIER GLASS AND DANIEL HAN-KWAN 



call hh cl unit vector, orthogonal to %. For / > 0, we denote 

-Hi ■.= n+[-l,l]nH- 
Since T-L is closed in T^, note that we can define d > such that 
(1.9) n2d C CO, 

and such that Ad is less than the distance between two successive lines in s~^{T-L). 

We first prove an "asymptotic" controllability result, for which the speed of light is consid- 
ered as a parameter. This next result will be valid for a large speed of light. 

Theorem 2 (With the strip assumption). Assume that U! satisfies the strip assumption. Let 
bo € //'^(T^;]R) be a magnetic field such that 

/ bodx^ 0. 

There exists Tq > such that for T > Tq, there exist k > and ct > satisfying the following. 
Ifc> CT, if ifo, Eo) G H^iT'^ X M2) X H^{T^) and fi € H^{T^ x M^) ^re such that fo, fi are 
compactly supported in v and satisfy 



(1.10) 


div Eo= fo dv 


- / fodvdx, 








(1.11) 


/ /o dv dx = 


/ fi dv dx, 




/t2xR2 


(1.12) 


II (/O, ^0)11/^3 <K, 


Il/l||/f3 < 1^, 


and if 






(1.13) 


Bo = 


cbo, 



then there exis ts a control G e H^{[0,T]xT^ xR"^), supported in uj, such that the corresponding 
solution {f,E,B) to the Vlasov- Maxwell system (|1.3p with initial data (fo, Eo, Bq) satisfies: 

f\t=T = /i- 

Remark 1.3. The scaling in (|1.13p is in fact natural if we want to recover the convergence 
of Vlasov-Maxwell to Vlasov-Poisson with an external magnetic field, in the classical limit 
c +oo: see Paragraph 12.41 

From the physical point of view, one can remark that Bo is in CGS units, namely Gauss 
(which is a theoretical physics unit), whereas — is in SI units, namely Teslas ("usual" unit). 
Thus, imposing that Bq = cbo is equivalent to imposing that the intensity keeps a non-zero 
value when c becomes large. So this assumption is also natural from this point of view. 

As a consequence of the previous result, we are also be able to prove the following theorem, 
in which the speed of light is considered as a fixed parameter, let us say c = 1. 

Theorem 3 (With the strip assumption). Let c = 1. Assume that uj satisfies the strip 
assumption. Let bo G i?'^(T^) be a magnetic field such that 

bo dx 0. 
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Then there exist Ti > and Ai > such that the following holds. For any A € (0, Ai], 
for any T > Ti/A, there exists k > such that, if (/o,-Eo) G H'-^iT^ x M^) x H'-^iJ^) and 
fi € H^(T'^ X M^) are compactly supported in v and satisfy 



(1.14) 


div Eo= fo dv 


- / fodvdx, 








(1.15) 


/ /o dv dx = 


/ /i dv dx, 








(1.16) 


||(/o,i?o)||//3 <K, 1 


|/l||//3 < K, 


and if 






(1.17) 


Bo{x) = 


A 60(2;), 



then there exists a control G S i?^([0, T] x xR^), supported in uj, such that the corresponding 
solution {f,E,B) to the Vlasov-Maxwell system (|1.3p with initial data {fo, Eq, Bq) satisfies: 

(1.18) f\t=T = h- 

Remark 1.4. The choice c = 1 in Theorem [3] is naturahy arbitrary (we could have taken any 
positive constant). Note also that taking A small corresponds to an assmnption of smallness 
on the initial magnetic field Bq. But this yields a result of controllability for a large time. As 
we will see in the proof, the effect of the magnetic field is indeed crucial to get the result. 

1.3. A short review on the Cauchy problem. Let us now provide a short review on the 
classical Cauchy theory for the Vlasov-Maxwell equation. To our knowledge, the very first 
results concerning these equations concern the local in time existence of classical solutions, 
and were obtained in independent works by K. Asano [1], P. Degond [9] and S. Wollmann 
|29| 130). All these contributions are based on Sobolev type function spaces. For a "sharper" 
theory based on spaces, we refer to the fundamental works of R. Glassey and W. Strauss 
|17| [T8l 120) . where the structure of the Vlasov-Maxwell equation (kinetic transport + wave) 
plays a crucial role. Concerning the global theory, there have also been several remarkable 
contributions: in 3D, for equations set on the whole space, R. Glassey and W. Strauss gave 
in the pionneering work |17| a famous criterion for global existence (see also the more recent 
contributions of S. Klainerman and G. Staffilani |21| . and F. Bouchut, F. Golse and C. Pallard 
[5] for alternative proofs); further studies around this criterion have been at the core of many 
papers (see for instance |19 | [25 1 [25]). but the unconditional global existence in three dimensions 
is still an open problem. In 2D, in the whole space case, the problem of global existence was 
settled by R. Glassey and J. Schaeffer |16) as well as for the so-called 2D- 1/2 case |15| . 

We also mention the fundamental contribution of R.J. DiPerna and P.-L. Lions, who built 
in |11| global weak solutions (in the energy space) to the Vlasov-Maxwell system, but we shall 
not use this type of solutions and will only consider classical solutions, for which characteristics 
are well-defined. 

1.4. Organization of the following of the paper. Let us describe how the paper is orga- 
nized. In Section [21 we explain the strategy of the proofs. In Sections [3] and [U we construct 
some relevant reference solutions, in order to apply the return method. Section [3] considers the 
case of Theorem [1] which corresponds to the situation with geometric control condition, while 
Section m deals with the case of Theorems [2] and [3] (with the strip assumption). For the latter, 
we will revisit the case of Vlasov-Poisson with external force fields of |14j and rely on the fact 
that the dynamics of Vlasov-Poisson with external magnetic field is a good approximation of 
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the dynamics of Vlasov-Maxwell, when the speed of hght is large. Finally, in Section [5l we 
present the approximation scheme which will allow to give a solution of the full non-linear 
problems and conclude the proofs of Theorems [U [2] and |3l This scheme involves an absorption 
procedure which yields some new technical difficulties (that do not appear when one considers 
the Cauchy problem for the relativistic Vlasov-Maxwell system without source). Finally we 
give in an Appendix the proof of a key lemma for Theorems [2] and [31 

2. Strategy of the proofs 

As a preliminary, we set up in this section the strategies that will be implemented in this 
paper. 



2.1. Some notations. Let us first give some useful notations. 

For T > 0, we denote Qt := [0,r] x x por a domain 0, we write C^(0), for / G N, 
for the set C\VL) n W^'°°{VL). For any / G N, we denote by H\T'^ x M?) the usual Sobolev 
space defined in x M^. 

For X in and r > 0, we denote by r) the open ball of center x and radius r, and by 
S{x, r) the corresponding circle. All radii will always be chosen small enough so that S{x, r) 
does not intersect itself. 

Let us also introduce some notations regarding characteristics and auto-induced electromag- 
netic field. Let F{t,x,v) some force field with Lispchitz regularity and sublinear at infinity in 
V. Let s>0, which corresponds to an "initial" time and {x,v) € x M^. 

We call {X{t,s,x,v),V{t,s,x,v)) the relativistic characteristics associated with F, the so- 
lutions of the system of ODEs: 



V, 



( dX 

~dt 

X{s, s,x,v) = X, V{s, s,x,v) = V, 



with V 



V 



1+^ 



We can also define {X{t, s, x, v),V{t, s, x, v)) the classical characteristics associated with F: 



dV 



dt 



X{s, s,x,v) = X, V{s, s, x,v) = V. 
In the following, when not otherwise specified, the characteristics under study are relativis- 



tic. 



Note that the characteristics (relativistic or classical) are well defined by the classical 
Cauchy-Lipschitz theorem. Usually, when there is no ambiguity, we will simply write (X, V) 
instead of the full notation {X{t,0, x,v),V {t,0, x,v)). 

Finally, for any distribution function f(t,x,v) sufficiently smooth and with sufficient de- 
cay in velocity, we will denote E^{t,x) and B^{t,x) the auto-induced force fields, which are 
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solutions of the Maxwell equations: 

dtE^ + ccurlS-^ = - / fvdv, dfB^ + ccurlS-^ = 0, 

divE^ = / fdv- / fdvdx, dw = 0, 
with some relevant initial data (depending on the context). 

2.2. The return method. Since we intend to show some controllability properties on the 
Vlasov-Maxwell system, for small data (that are close to 0), it would be natural to apply the 
following proof scheme: 

(1) Linearize the equations around and prove some controllability properties on these 
linear equations. 

(2) Prove that these properties are somehow preserved by the full equations, by using 
some fixed point theorem. 

Unfortunately, we are in a situation where this scheme fails. If we consider the linearized 
Vlasov-Maxwell equation, around the trivial state (/, E, B) = (0, 0, 0), we obtain the following 
equation: 

(2.1) dtf + v-V^f = 1^G, 

which is the relativistic free transport equation. The characteristics associated with the trans- 
port part are simply straight lines. By Duhamel's formula, we obtain the explicit representa- 
tion for /: 

f{t,x,v) = fo{x-tv,v)+ / {l^G){s,x - {t - s)v,v)ds. 

Jo 

One can observe that there are two types of obstruction to controllability for the (linear) 
relativistic free transport equation. 

• (Small velocities) The first obstruction concerns the small velocities. If the velocity v 
of a charged particle driven by (j2.ip is not large enough, then it will not be able to 
reach the control zone u in the desired time T, so that {l^G){s,x — tv,v) = for all 
t € [0, r]. In particular, if c is not large enough with respect to the desired control 
time, then this obstruction appears for many trajectories. 

• (Bad directions) The second obstruction is of geometric type: it is possible that a 
particle has initially a bad direction, so that, even if the modulus of its velocity and 
the speed of light are very large, it will never reach the control zone. Thus we cannot 
influence the value of the distribution function along this trajectory. There is no such 
problem when u satisfies the geometric control condition, but this problem can occur 
under the strip assumption. 

To overcome this difficulty, we shall rely on the return method of J.-M. Coron (see the book 
[8] for many examples of applications of this method for several non-linear PDEs, especially 
coming from fluid mechanics). The idea is to build some relevant reference solution of the 
full PDEs (with a well-chosen source supported in the control set w), starting from the trivial 
state and reaching again after some time fixed in advance, around which the linearized 
equations do enjoy nice controllability properties (as opposed to the trivial solution 0). Then, 
in the fixed point scheme designed to construct a solution of the non-linear problem, we will 
linearize the PDEs around this reference solution, and not the trivial one. 
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Figure 1. The two obstructions for controllabiUty 

For our problem, the key argument is to find a reference solution (/, E, B) to the relativistic 
Vlasov-Maxwell system, with some suitable source in w, starting from (0, 0, 0) and coming back 
to (0,0,0), and such that the characteristics associated with E + satisfy: 

(2.2) G T2,Vt; € IR^3^ G [0,r], X{t,x,v) u. 
Note indeed that the linearized equation reads: 

- v^- 

dtf + v- Vxf + {E + —B) ■ V,f = l^G, 

c 

so that (j2.2p means that the two above obstructions do not occur. 

It is clear that this condition on the characteristics is necessary to get some controllability 
properties on the linearized equations; we will show later that this will also be sufficient to 
get local controllability for the full equations. 

Actually, in view of the fixed point procedure which we use to build a solution of the 
non-linear problem, we will need to prove a slightly refined version of ()2.2p : 

(1) In the case with GCC, we consider open sets uj' ,oj" contained in cj, which still satisfy 
GCC, such that 

w' C a7 C C a/' C w. 
We ask that the characteristics associated with E + ^B satisfy: 

(2.3) Vx G t2,Vu G M^3^ G [0,T], X{t,x,v) ^ oj' . 

(2) In the case with the strip assumption, we consider G a; and tq > such that 
B{xq, 2ro) is contained in w. We ask that the characteristics associated with E + —B 
satisfy: 

(2.4) Vx G T^V'^; G M^3^ G [0,r], X(t, x, t>) G S(xo, ro). 

2.3. The case of Theorem [H With geometric control condition on the control set we 
can rely on exact controllability properties for the Maxwell equations, which were already 
studied by K.-D. Phung in |24| (actually in the more involved context of open sets of with 
boundary) . This controllability tool will be given in Theorem [H The crucial point is that 
there exist stationary solutions of Maxwell which are relevant for (j2.3p . This will allow us to 
build a suitable reference solution as follows. Starting from zero, by using a suitable control 
G, one can emit a wave from so that the electromagnetic field reaches such a stationary 
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solution. Then one stays at this point for some time interval, and then returns to zero by 
emiting again a wave from u. 

2.4. The case of Theorem [2j, Without the geometric control condition, it would be at 
first sight tempting to try to use some results based on the approximate controllability of 
wave equations (such results are for instance given for the wave equations in the survey of 
D. Russell |26| and in the book of J.-L. Lions |22|). Unfortunately, there are at least two 
obstructions. First, it is not clear at all that such results hold for the Maxwell equations 
(because of some consistency issues, see equation (|3.ip below, that do not appear for the 
wave equations). Furthermore, such a strategy is likely to fail, since one has to find a way to 
"sustain" for some time the appropriate electromagnetic field {E, B), so that the trajectories of 
the particles get sufficiently affected by them. This was possible with the geometric condition 
condition on u, only because it is possible to exactly reach some stationary solution of the 
Maxwell equations. With only approximate controllability, it seems no longer possible to do 
so. 

Instead, we propose a strategy based on the large speed of light asymptotics (in other words 
the non-relativistic limit) of the Vlasov-Maxwell equation. From the physical point of view, 
this limit is well understood: it is well-known that as the speed of light goes to infinity, the 
Maxwell equations are approximated by the Poisson equation: this corresponds to the so-called 
electrostatic approximation, which is valid when the speed of light is large enough compared 
to the typical velocities of the system. Concerning the transport part, that is, the relativistic 
Vlasov equation, one can readily see that the effect of a large speed of light is to lead, at least 
formally, to the "classical" transport setting, that is: 

V 

V. 



Y _(_ 1"!^ c-i.+oo 

When c -|-oo, one can indeed rigorously prove the convergence of classical solutions of the 
Vlasov-Maxwell system towards classical solutions of the Vlasov-Poisson system. This was 
done independently in works by K. Asano and S. Ukai [2], P. Degong [9] and J. Schaeffer |27| . 
in the framework of equations set in the whole space M^. We recall that the Vlasov-Poisson 
(without magnetic field) system reads: 

(dtf + v- VJ + div„ [V,.(/./] = 0, 



Ax(p = J f dv — J fdvdx, 
, f\t=o = fo- 

Somehow motivated by the convergence of Vlasov-Maxwell to Vlasov-Poisson, we intend 
to show that the study of the controllability of the Vlasov-Poisson system with external force 
fields, which was done in our previous work |14| . can be useful in the Maxwell context. This 
means that we rely on the fact that the Vlasov-Maxwell dynamics is close to that of the simpler 
Vlasov-Poisson with external magnetic field system, when c is very large. The crucial points 
are Lemma 14.21 (which shows how magnetic fields have a "bending" effect on characteristics 
and help to overcome the bad directions) and the Approximation Lemma 14.31 which allows 
us in some sense to study a Poisson equation instead of Maxwell equations. To apply the 
approximation lemma, some difficulties appear, due to some consistency issues in the Maxwell 
and Poisson equations, and we will rely on the strip assumption on the control set to to solve 
them. 
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Remark 2.1. Taking Bq = cbo in ()1.13p . where c is the speed of light, is the "right" scaling 
if one wants to recover the Vlasov-Poisson system with an external magnetic field bo in the 
classical limit c — )■ +00 (taking Bq = Ijq yields no magnetic field in [2i 



Remark 2.2. In some sense, the spirit of our proof is quite similar to that of the proof of 
approximate controllability of the incompressible Navier-Stokes equations with Navier slip 
boundary conditions, by J.-M. Coron [7], whose basic principle relies on the fact that for small 
viscosities, solutions of the incompressible Navier-Stokes system behave like solutions of the 
incompressible Euler equations. 

2.5. The case of Theorem [Sj scahng properties of the Vlasov-Maxwell system. 

Finally, for the proof of Theorem [3l we will also exploit some scaling properties of Vlasov- 
Maxwell, that we recall here. One can observe that the Vlasov-Maxwell system is "invariant" 
by some change of scales. This is a kind of generic property to Vlasov equations. More 



precisely, when (/, E, B) is a solution of Vlasov-Maxwell in [0, T] 
(2.5) f^{t,x,v):= f{\t,x,v/\), 



X R^, then for A / 0, 



(2.6) 



E^{t,x) := X^E{Xt,x), 



(2.7) B^{t,x) := X^B{Xt,x). 

is still a solution in [0,T/A] x x of a rescaled Vlasov-Maxwell system: 



(2.e 



dtr + V ■ V,./^ + div. 



dtE^ + cXcuilB^ 



1 



{E^ + —vAB^)f 
cX 



f^vdv, 



dtB^ + cX curl E^ 



0, 



div E' 
/|t=o 



A 



f^dvdx, div 5^ = 0, 



E\ 



t=o 



E, 



A 
5 



Bi 



t=o 



B, 



' 



with V 



1+ 



b|2 



This implies that choosing a large parameter A S> 1 artificially increases the speed of light. 
This is the key point to deduce Theorem [3] from Theorem [2j Another important choice is 
A = — 1, which shows that the equations are reversible. This will allow to simplify the proof 
of Theorem [H we will only prove the result for /i = outside uj, {Ei, Bi) = (0, 0) and deduce 
the general case by reversibility (actually it is possible to avoid this argument, see Subsection 

In another direction, we mention that the choice < A ^ 1 can be useful to prove global 
controllability results for Vlasov systems (that correspond to results without restriction on 
the size of the data) , and was indeed a crucial argument in the Vlasov-Poisson case in \13\ I14j . 
Unfortunately, we will not be able to do so for this problem, since it is not compatible with 
our method using the bending effect of the magnetic field. 

3. The reference solution with geometric control condition 

In this section, we assume that the geometric control condition is satisfied by the control 
set UJ, and we build the reference solution (f,E,B) which is central in the proof of Theorem 
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[TJ To this purpose, we first describe in Subsection 13.11 some results on the controUabihty 
of the MaxweU equations. These results are needed for the construction of the reference 
solution, which is given in Subsection 13.21 (in particular we will explain how to obtain a 
relevant distribution function /). 

3.1. Controllability of the Maxwell equations with the geometric control condition. 

The controllability problem which we consider in this paragraph is the following one. Let T > 
be some control time and let {Eq, Bq, Ei, Bi) be some electromagnetic field (belonging to some 
appropriate function set), the question is to find j and p supported in co such that 



f dtB + ccurl^ = 0, dtE + cmrlB 



-J, 



divE 

E\t=o ~- 



-- P- 

--El, 



pdx, div B = 0, 



J2 

B\t=o = 

B\t=T 



Bo, 
--Bi. 



Using the reversibility of the Maxwell equations, it is straightforward to see that this is equiv- 
alent to the next controllability problem. Let {Eq,Bq) be some electromagnetic field, find j 
and p supported in lo such that 



' dtB + ccurlE' 
div^ 
E 



El 



^\t=o -- 

\t=T 



-- P- 

Eq, 
-0, 



B\t=T 



0, dfE + ccmlB 
pdx, div B = 0, 

■ Bq, 

0. 



-J, 



This controllability problem was solved by K.-D. Phung in |24| for equations set in a smooth 
domain of . One difficult point comes from a consistency issue. To solve the Maxwell 
equations, the source functions have to satisfy the local conservation of charge: 



(3.1) 



9tp + V^. • j = 0, VtG [o,r], VxGT2. 



K.-D. Phung bypasses this difficulty by showing that it is possible to find a control such that 
p = and divj = and we follow this approach. We will actually use a slight variation on 
Phung's theorem. Although we will not give a complete proof of this result, we provide for 
the sake of completeness several lemmas, which, when put together, yield the claimed result. 

Using the celebrated H.U.M. (Hilbert Uniqueness Method) of J.-L. Lions, which consists of a 
duality argument, see |22| . proving the controllability is equivalent to proving an observability 
inequality: 

Lemma 3.1. The Maxwell equations can he controlled in time T if and only if the observability 
inequality on the adjoint system is satisfied: there exists C > such that for any Eq,Bq € 
L^(T^), with V • -Eq = 0, the solution {E,B) to the homogeneous Maxwell system 



(3.2) 



{ dtB + ccuvlE = 0, 
div£; = 0, divB 



El 



\t=o 



En 



B 



\t=o 



dtE + ccMilB 
0, 

= -Bo, 



0, 
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satisfies 



\E\ 



L2([0,T],L2(T2)) 



+ ll^l 



L2([o,T],L2(T2)) 



< c 



\E\ dxds. 



Now we state another lemma due to Phung. This result actually comes from the conserva- 
tion of the electromagnetic energy and a decomposition of the electromagnetic field using a 
suitable vector potential. 



Lemma 3.2. There exists Ci > such that for any Eq,Bq € L^(T^) 
solution {E, B) to ()3.2p satisfies: 



with V • £'o = 0, the 



I^IIl2([0,T],L2(T2)) < C'i|I^IIl2([o,T],L2(T2))- 



In particular, it is sufficient to prove observability for the electric field only. 

Then, writing the equation on E as a wave equation, K.-D. Phung proved that the ob- 
servability inequality is true in a domain with boundaries, using results of propagation of 
singularities due to L. Hormander, R. Melrose and J. Sjostrand, as in the papers of J. Rauch 
and M. Taylor [25] and C. Bardos, G. Lebeau and J. Rauch 1^. 

Lemma 3.3. The observability inequality holds true ifT is large enough. 

In view of the control of the Vlasov-Maxwell system, we additionally need to build smoother 
controls, so that the associated electromagnetic field (E, B) remains smooth (at least with 
Lipschitz regularity) for all times. This can be done by using the method of B. Dehman and 
G. Lebeau |10| (see also the paper of C. Bardos, G. Lebeau and J. Rauch [5]) or by applying 
the systematic procedure proposed in the paper by S. Ervedoza and E. Zuazua |12| . Gathering 
all results together, we obtain the following theorem. 

Theorem 4. Assume that the control time T > is large enough. Let /c € N*. For any 
Eo,Bo,Ei,Bi G H''{T^), with div Eq = div^i = 0, there exists a control function j € 
n';^QH'{[0,T];H^^'{T'^)), supported in [0,T] X oo, satisfying j\t=o = j\t=T = 0; such that for 
all t G [0,T], div J = and such that the solution {E,B) to the system: 

dtB + ccmlE = 0, dtE + ccuilB = —j, 

div^ = 0, divB = 0, 



satisfies E\ 



t=T 



E, 



Bi 



E\t=o — Eq, 
t=T = Bi. 



B 



\t=0 



Bn 



3.2. Construction of relevant reference solutions. Let uj', uj" be some smooth open sets, 
contained in w, which still satisfy GCC, such that uj' C uj" C uj" C uj. Our aim is now to apply 
the controllability results for the Maxwell equation, in order to build, for some sufficiently large 
T > 0, a suitable smooth reference solution (/, E, B) for the full Vlasov-Maxwell equations: 



dtf + v- Vccf + div„ 



{E + —B)f 
c 



dtE + c curls 



fvdv, dtB + ccuilE = 0, 



divE 



fdv 



R2xT2 



/ dv dx, 
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where G is a suitable source and which satisfies the following properties: 



(3.3) 
(3.4) 



/|i=o = 0, ^|t=o = 0, 5|<=o = and f^t^^ = 0, ^|t=T = 0, B^t=T = 0. 
For all t G [0,T],v G M^7(^, . , v) is compactly supported in oj. 

- v^- 

The characteristics [X, V) associated with E -\ B satisfy the property: 



(3.5) 



for any x G T^u G there is t G [T/9,8r/9] such that: 

X{t,0,x,v) G u)', and \V\{t,0,x,v) > 4. 

The main result of this section is the following: 

Proposition 3.1. //T > is large enough with respect to c, there exists a reference solution 
{7,E,B) such that 1^, 1^, 1^ hold. 

Proof of Proposition \3.1\ To build such a reference solution, as explained in Section [21 we have 
to overcome the obstruction coming from slow particles. The idea is to use an electric field to 
accelerate particles whose initial velocity is too low. 

The elementary (yet crucial) observation is that on the torus, any constant electromagnetic 
field {Ei,Bi) G X M is a stationary solution of the Maxwell equations (without source). 

1. In a first time interval [0, Ti] C [0,T), we take as the reference solution 



which is of course a trivial stationary solution of Vlasov-Maxwell; by compactness, thanks to 
the geometric control condition satisfied by uj' , there exists Ti > such that for all a; G T*^, 
and ah u G such that \v\ > c/2, 



Then there exists m > 4, such that for any \v\ > m, \v\ > c/2, which allows to say that ()3.5p 
is satisfied for all particles whose initial velocity v has a modulus larger than m. 

2. There remains to take care of all other particles, whose velocity is not high enough. 
For that, the idea is to use the control for the Maxwell equations (Theorem [4]) to steer the 
electromagnetic field from (0,0) to (Si,0), with Ei := (1,0), to wait for a sufficiently large 
time so that all particles get a sufficiently high velocity and meet the control zone as in the 
first step, and finally, to use once more the control for Maxwell to bring the electromagnetic 
field back to (0,0). Let us now make this construction explicit. 

Let T2 > be a large enough time such that we can apply Theorem |4] and find a smooth j2 
(with V • j2 = 0), supported in u such that the following holds: 



{f,E,B) = {0,0,0), 



(3.6) 



X + tv G uj' , for some t G [0, Ti]. 



( dtB + ccurlE = 0, dtE + ccuA B = -j2, 



div ^ = 0, 



dwB = 



< 



El 



t=Ti 



= 




< E\t=Ti+T2 
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Now, we can consider two functions Zi G Cq 
constraints: 

> in M^, 



(for i = 1 or 2) satisfying the following 



(3.7) 



We then set 



Supp Zi C 5(0, 1) 
Zi(v) dv = 0, 



Zi{v)vdv = {5i=i,5i=2)- 



f = Zi {v)jl + Z2 {v)j2 , where j2 = (ia ,^2)1 



on [Ti,Ti + T2]. We take {E,B) = (E,B) and we notice that {f,E,B) is a solution of the 
Vlasov-Maxwell system on [ri,Ti + T2], with a suitable source (supported in uj): 



G{t, X, v) := dtf + v- +[E + —Bj ■ V,/. 

Notice that G is compactly supported in uj, by construction of j2 and /. 

The effect of constant the force field Ei is to accelerate particles (even if it takes a very long 
time) : denoting by {X, V) the characteristics associated with Ei , we can prove the existence 
of some T3 > such that for any (x, z;) € x M^, there exists t G [Ti + T2, Ti + T2 + T3] , such 
that: 



(3.8) Vx G T^ \/v £ s.t. \v\ < m, 



X{t,0,x,v) e uj' and \V\{t,0,x,v) > a. 



Then, on the interval of time [Ti + T2,Ti + T2 + T^], we consider (/, E, B) = (0, 0) as 
the reference function, which is a stationary solution of Vlasov-Maxwell without source. 



3. Finally, one finds a smooth J3 such that the following holds: 

+ ccurl^ = 0, dtE + ccm\B = 
div^ = 0, divi? = 0. 



-J2, 



i=Ti+T2+T3 - 
t=Ti+2T2+T3 



El, 
= 0, 



B\ 



t=Ti+T2+T3 



B 



\t=Ti+2T2+T3 



0, 

: 0. 



and we can then define 



f = Zi {v)jl + Z2 {v)ji , where J3 = {jl , j| ) , 

on [Ti + T2 + T3,Ti + 2T2 + T3]. As before, we take {E,B) = {E,B). Finally we choose 
T > Ti +2T2 + r3. By construction, (|3.3p and ()3.4p are satisfied. We observe finally that (j3.5p 
is satisfied for some times in [0,T], because of p.6p and (|3.8p . Obtaining the time interval 
[T/9, 8T/9] in (|3.5p is a matter of translating the solution in time and extending it by when 
necessary. 

□ 



4. The reference solution under the strip assumption 

In this section, we build a reference solution (/, E, B) for the Vlasov-Maxwell system, in 
the case where we make the strip assumption on the control set. We shall also assume that 
the speed of light satisfies c > cq (where cq > will be fixed later, see Lemma l4.2p . This 
construction is central in the proofs of Theorems [2] and [3l 
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This section is divided into three sub-sections. First, we build a relevant reference solution 
for a relativistic Vlasov-Poisson system. Then, we state an approximation lemma for the 
Maxwell equations (as c — > -|-oo), which we use in the last subsection to build a reference 
solution for the Vlasov-Maxwell system. 

4.1. The case of Vlasov-Poisson v^ith an external magnetic field. Let xq £ uj and 

ro > small enough such that B{xq, 2ro) C OJ. Let b some nonzero constant external magnetic 
field. Let R he a large enough positive constant; we construct a reference solution /, which 
depends on R, and which will be used to solve the controllability problem for data /o and fi 
supported in x B{0,R/2). 

We first aim at finding a reference solution / for the relativistic Vlasov-Poisson equation, 
with the external magnetic field b: 



(4.1) 



dtf + v- V,./ + div„ [(V,(^ + v^b)f 
A^Tp = p- 1, 



where G is a suitable source and which satisfies the following properties, for some large enough 
T > 0: 

(4.2) 7|i=o = l\t=T = 0. 

(4.3) Vt G [0,r],Vt; G 7(t,. , v) is compactly supported in u. 

(4.4) The characteristics (V, V) associated with \7Tp + v'^b satisfy the property: 
for any x G T^u G B{0,R), there is t G [r/9,8T/9] such that 

X{t,0,x,v) G B{xo,ro/2), and \V\{t,0,x,v) > 5. 

Building such a reference solution has been already done in |14| for the classical Vlasov- 
Poisson equation; actually it was done for any arbitrarily small control time T > 0, but the 
generalization to this relativistic version of Vlasov-Poisson imposes to take a large enough 
control time. Furthermore, having in mind an application to the Vlasov-Maxwell equation, 
because of consistency issues, we have to impose two additional conditions. First that the 
local conservation of charge is satisfied everywhere, 

(4.5) dtp + divj = 0, Vx G T^, Vt G [0,T], 

(where p := J f dv and j = J fvdv) and second that the "zero-mean current" condition holds 

(4.6) / jdx = 0, yte [0,T]. 

Let us introduce some notations before getting to the construction of /. 

Notation. We denote hy T> a line of T^, which does not cut the control zone uj (reduce 
UJ if necessary, but in a way that it still satisfies the strip assumption) and n a unit vector, 
orthogonal to V. 

We also recall that ?^2(i is defined in (|1.9p . 



We have the following lemmas, which are taken from [T3], with modifications to ensure that 
()4.5p and (|4.6p are satisfied. The first lemma states that we can find some relevant electric 
field, whose effect will be to accelerate particles. 
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Lemma 4.1. Let b G M 6e given. Let r > and M > 0. There exists M > 0, £ e 
C°°([0,r] X T2;M2) and if G C°°([0,r] x T'^;R) which do not depend on c, satisfying 

(4.7) £ = -Vip tn[0,T]x{T^\nd), 

(4.8) Supp{£) C {0,t) xT'^, 

(4.9) A^p = m [0, r] x (T^ \ Hd), 

(4.10) / Vip-ndx = 0, 

Jv 

such that, for any ^ G L°°{0,T;W'^'°°{T'^ x M^)) satisfying \\d\\L°- < 1, «/ (^, V') are i/ie 
characteristics corresponding the force ^ + £ + ij^b, then 

(4.11) y{x,v) X B{0,M), V{t,0,x,v) £ B{0,M)\B{0,M + 1). 

The proof is close to that of \14:\ Proposition 5.2]. The difference with the statement of |14) 
appears in (|4.10p . which will be crucial to ensure (|4.5p . To prove all requirements, one has to 
find w G C°°(T2;]R2) satisfying : 

(4.12) curl-u; = in 

(4.13) divw = inj'^ \ nd, 

(4.14) \w{x)\ > for any x in \ Ud, 



(4.15) / w-ndx= / w ■ dr = 0. 

Jv Jv 

Such a function was constructed in |13| Appendix A. 2]. Once it is obtained, one can take 
£ := w and observe that it coincides with a gradient outside of Tid- 

Note in particular that the solution here does not correspond to the one giving rise to the 
force field n (which is the "most natural" one could think of; it satisfies (j4.12p . (|4.13p and 
()4.14p but not (j4.15p ). The effect of this force field £ is only to accelerate particles, and it is 
very unlikely that it allows to "cure" the bad directions. 

For that issue, we will fully rely on a "bending" effect of the magnetic field b, which is 
described in the next lemma. 

Lemma 4.2. Let M > 1. Given b ^ 0, there exist: 

• Co > depending on b,Xo,?'0; 

• 777- > depending only on b,XQ,rQ, 

• Tq > depending on b, XQ,rQ and M, and 

• K depending on b, XQ,rQ and M, 

such that for all J G L°°(0, T; Tyi'°°(T2 x M^)) satisfying < k, if c > cq then the 

characteristics {X,V) associated with ?)^b + 5^ satisfy: 

(4.16) 

Vx G T^Vu G such thatM> \v\ >m,3te (ro/4, 3ro/4), X{t,0,x,v) G B{xo,ro/2) 

and for all s G [0, Tq], -Li < \V{s, 0, x,v)\ < 2\v\. 



The proof of Lemma 14.21 (actually, of a generalized version of Lemma 14. 2p is given in the 
appendix. 
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With these two lemmas in hand, we now have: 

Proposition 4.1. There exists T > large enough, such that there exists a reference solution 
(7,^) such that ()42]) . (fO]l . (|4ll|) . (fi3]l . ()46]) are satisfied. In addition, p := J f dv and 
j := J fvdv do not depend on c. 

Proof of Proposition \4-l\ The proof is divided into three steps. 



1. First, we define the reference potential : [0,T] x ^> M as follows. We apply Lemma 
14.21 for M = R+1, and we obtain some time Tq and some m > 10 such that (|4.16p is satisfied. 
Then we apply Lemma [4.11 with r = Tq and 

(4.17) M = max(^m + l,32ro(|b| + l)), 

and obtain some (^27 ^2 and some M > such that ()4.1ip is satisfied. We apply Lemma 
for M = M , and we obtain some time Ti (with the same m > 0). We set 

OfortG [0,ro]U[2ro,2To+Ti], 
^2(i-To,-) for t G [ro,2ro]. 



and 

We define T := 2To + ri. 



Ofor t € [o,To]u[2ro,2ro+ri], 

£2{t-To,-) forte [ro,2To]. 



(4.18) 



2. Let us now introduce a first distribution function /. Consider a function Z € C^(M^;M) 
satisfying the following constraints: 

r Z > in 

Supp Z C Sr2(0, 1), 

/ Z{v)dv = l, / Z{v)vdv = 0. 
We introduce / = f{t, x, v) as 

(4.19) fit,x,v):=Z{v)ATp{t,x). 

Of course, / satisfies ()4.1I1 in [0,T] x x M^, with source term 

(4.20) G{t, X, v) := dj + v ■ Vj + (v^ + bv^^ ■ V J, 

which is supported in [0, T] x B{xQ,rQ) x M?. Up to an additive function of t, the function 7p 
satisfies the Poisson equation corresponding to /. We can observe that / also satisfies: 

f\t=o = 0, f]t=T = 0. 

Finally, we denote 



(4.21) p{t,x):= / f{t,x,v)dv = Alf{t,x). 

By construction, according to Lemmas 14. II and 14.21 ()4.2p . (|4.3p and (j4.4p hold. 

3. Consequently, up to now, the solution (/,^) satisfies all requirements but (|4.5p and (|4.6p . 
During the interval of time [0,Tq] U [2Tq,2Tq + Ti], these are trivially satisfied since / = 0. 
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The problem appears only for the interval of time [To,2ro], where in general, we do not have 
the compatibility conditions: 



dt [ fdv + w^- [ 

JR2 Jm? 



fv dv = 0. 



and 

/ / fv dv dx = 0. 

To ensure these conservation laws, the idea is add a correction to /. Let 

g := ui{t,x)Zi{v) + U2{t,x)Z2iv), 

with u{t,-) := {ui,U2) compactly supported in an open set U such that U C OJ, where we 
recall that Zi (for i = 1,2) is defined by: 

r > in 

Supp Z, C B{0,1), 

Zi{v)dv = 0, / Zi{v)vdv = {5i=i , 6i=2 ) ■ 



Let us show that we can choose u such that: 

dt / fdv = -Vx ■ / gvdv, 

and 

/ g V dv dx = 0. 

For simplicity, we denote h := —dt J fdv. The problem is equivalent to finding a vector 
field u G C°°(M"'" x T^), satisfying all four conditions: 

(4.22) u\t=o = u\t=T = 0, 

(4.23) divu = h, Vt € [0, T], x G T^, 

(4.24) Vt G [0,r],Suppti(t, •) C w, 

/ udx = 0, VtG[0,r]. 
Jt2 

To build such a function u, we introduce 6 G C°°([0, T] x T^) a solution of the elliptic equation: 

A9 = h, for t G [0,r], x G T^. 

Observe that fj2 hdx = due to (|4.2ip . With the notations of Lemma 14.11 we notice that we 
can take 

e = -dtif. 

We thus conclude that V0 is divergence free outside u. Therefore, according to the Hodge- 
Poincare lemma (recall that uj contains a hyperplane of the torus with normal vector n), there 
exists ^' G C°°([0,r] X T^) and Q G M such that: 

V0 = V-^^' + an in \ uj, 

where a = V9 ■ ndx. But the reference potential ip defined above is designed as to satisfy 
Jjj'Vf - ndx = (see ()4.10p in Lemma l4.ip . This yields that a = 0. Now, we introduce a 



(4.25) 
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smooth cut-off function rj such that r] = on U, and r] = 1 on \ uj. We define a vector field 
u as: 

u = Ve- V^(r?^'), 

and one can finally check that u satisfies the four above conditions ()4.22p . (|4.23p . (|4.24p and 
(225]). 

We can now set: 

l = f + 9- 

Note that by construction, we have: 

Jdv = / fdv, 

SO that the electric field created by / (through the Poisson equation) is the same as the one 
created by /, and that 

J f vdv = J 'gv dv. 

One can readily check that / is a solution of the Vlasov- Poisson system with a suitable source 
in u, and that (|4.2p . (j4.3p . (j4.4p . (|4.5p . (j4.6p are satisfied (up to shifting and rescaling the time 
interval [0,r]). The fact that Jj and j do not depend on c (but only on cq) is due to the fact 
that the constructions of Lemmas 14.11 and 14.21 do not depend on c (but only on cq). □ 

4.2. The approximation lemma. Now, we want to use the following approximation lemma, 
which is valid for large values of the speed of light. The goal is to approximate the Vlasov- 
Maxwell dynamics by the simpler one of Vlasov-Poisson with an external magnetic field. 

For any function j4 : T'^ — ?> M, c? = 2 or 3, at least in L'^, we denote by the k-th Fourier 
coefficient of A, for k ^ Z'^. 

Lemma 4.3. Let d = 2 or 3. Let Eq, Bq, j, p some C°° functions. Let {E,B) the solution to 
the Maxwell equations: 

' dtB + c curl ^ = 0, dtE + c curl B = -j in [0, T]xT'^, 

divE = p-jpdx, divB = in[0,T]xT'^, 

{Eit=o = Eo, B\t=o = Bo znT-^. 
Assume that at initial time the compatibility conditions are satisfied: 

divEo = yo(0) - j piO) dx m [0, T] x T'^, 

divSo = inT'^. 
Assume that the local conservation of charge is satisfied for all times: 

Vx € dtp + V, • j = 0, 

as well as the zero-mean current property: 

(4.26) / jdx = 0. 

Jt2 
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Let be the solution of the Poisson equation: 

' curlE^oo = 2n [0,r] x T"^, 
div^Joo = yo - / pdx in [0,T] x T"^. 



Then, we have for all t > 0: 



(4.27) 

and 
(4.28) 



B- Bndx-B 



< 



C,,j{t + 1) 



L°°([0,t]xT<^) 



Eodx-E 



< 



c;,,(i+i) 



L°°{[0,t]xTd) 



where Cpj and C'^ j are constants depending only on p and j and where B and E are defined 
by their space Fourier coefficients: 



(4.29) 



B =0, 

_ kAE''{0) 

E =0, 



sm{tc\k\) + ^^(0) cos(te|A;|) for /c / 0, 



^ E = _ ^'AB'j(o) sin(tc|fc|) + [^^=(0) - 3ip''(0)\ cos(te|A:|) for k ^ 0. 



Remark 4.1. The initial data in this lemma are somehow ill-prepared. In contrast, in the 
well-prepared case, which would correspond to the additional assumptions: 

curli^o = 0, curli?o = 0, 

then we would have E = and B = and this would yield strong convergence results in 
L°°([0, t] X T''). In the present case, the convergence is only weak in time. 

Proof of Lemma \4-3\ The electromagnetic field E and B satisfies the wave equations: 

r d^E - c'A^E = -c^Vp - dtj, 



d'^B-c^A^B = ccurlj. 



Thus we have: 



r d^E^ + c'lkl^E" = -ikc^p^ - dtj", 
[dfB'' + c^\k\'^B" = cikAj\ 

that we easily solve, for k ^ 0, with Duhamel's formula: 

(4.30) 



E'' = —dtE^iO) sin(te|A:|) + ^^=(0) cos 



c\k 



{tc\k\)+ [ 
Jo 



c^ikp^s) - dtj^] 



..fc sin((t - s)c|A;|) 



c\k\ 



(4.31) B^ = —dtB^{Q)s\n{tc\k\) + B^{<d)cos{tc\k\)+ cikAf{sj 

c\k\ Jq c\k\ 



■fc/ Asm((t - s)c\k\) 



ds. 



ds. 
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Let us first deal with E . We have, by integration by parts: 



d'ikp {s) 



We remark that: 



sin((t-s)c|fc|) 



c\k\ 
—ik 



ds 



\k\ 



ik ^i., 



* ik 

1^' 



^p^{t) + j^p^{Q)cos{tc\k\)+ I j^dsp''{s)cos{{t- s)c\k\)ds. 



—ik 

w 

By another integration by parts, we get: 
/ -—:TdsP^{s)cos{{t-s)c\k\)ds 

Jo \kr 

ik 



p\t) = EUt)- 



c\k{' 



ik 



dtp^i^) sin(tc|fe|) + / ^r^dssP^is) sin((t - s)c\k\) ds. 



Moreover, we have: 



ik 

c|fcp 



dssp'^is) sm{{t - s)c\k\) ds 



c\k\ 



1 t 



Finally, in a straightforward manner, we have 



c\k\ 



< 



so that for k 0: 
(4.32) E''{t)-Ei{t) 



1^1 



m?\\Lt 

c\k\ 



ik 



sin(te|A;|) + ^'=(0) - tttoP i^) cos(tc|A:|) 



< 



T^\\dttp'\\Lr + 



\k\' 



1 WWIl^ Wfhr ] mj^L? 

C [ |fc|2 + \k\ ] ^ C\k\ 



For k = 0, one has, using the zero-mean current condition ()4.26p . that 

E\t) - Ei^it) = E\0). 

Finally we have: 



E-Eoo- I Eodx-E 



L°°([0,i]xT2) 



^tp^h 



t 



\k\ 



\k\ 



2 +nirJ+R"^*^''"^*°^ 



The summability is a simple consequence of the smoothness of the sources p and j. Hence, 
the estimation (|4.28p holds with 



1 / mp^'u? 

c 



\k\ 



1 

W\ 
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The magnetic field B is treated likewise. 



□ 



4.3. The reference solution. We now rely on the approximation lemma to find a suitable 
reference solution for the Vlasov-Maxwell system, provided that the speed of light is large 
enough, and that the initial magnetic field Bq is of the form c bo where bo satisfies 

60 dx 7^ 0. 

More precisely, we look for c large enough such that there exists a smooth reference solution 
(/, E, B) to the Vlasov-Maxwell equation: 



(4.33) 



dtf + v ■ V^f + div^ 



dtEi - cdxnB 



{E + —B)f 

c 



fvidv, dtE2 + cdx^B 



fv2 dv, 



dtB + d,,E2 - d,,Ei = 0, 



fdv 



f dv dx, 



T2xl 



with G is a suitable source, which satisfies the following properties: 

(4.34) 7\t=o = 0, J\t=T = 0, E\t=o = and B\t=o = Bq, 

(4.35) yt € [0,r], G R^, J{t,. , v) is compactly supported in to, 

- v^- 

(4.36) The characteristics {X, V) associated with E -\ B satisfy the property: 

for any x € T^u G there is t G [r/10,9r/10] such that 

X{t, 0, X, v) G B{xo, ro) and \V\(t, 0, x, v) > 4. 

Proposition 4.2. Let T > large enough. There exists c{T) > such that for any c > c{T), 
there exists a reference solution {f,E,B) of (j4.33p such that (|4.34p . (j4.35p and (j4.36p hold. 

Proof of Proposition \4.2\ Let (/, V^) the reference function in Proposition 14. 11 corresponding 
to the case of Vlasov-Poisson with external magnetic field, relevant for some large enough 
T > 0. Let {E, B) the electromagnetic field satisfying Maxwell equations with p = J f dv and 
j = J fvdv as sources; we recall that by construction, they do not depend on c, but only 
on the lower bound cq (this fact is crucial). We observe that {f,E,B) is a solution of the 
Vlasov-Maxwell system (with a suitable source in u). In addition, (|4.34p and (|4.35p clearly 
hold. 

Let us define (recall (jl.l3p ) 

b = - I Bodx = I bo dx. 

We introduce three different characteristics: 

• the characteristics {X,V) are associated with \/lp + v-^b, 

• the characteristics {X, V) are associated with E — E + ^{B — B), 

• the characteristics {X,V) associated with E +—B. 
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Let now show that if c is large enough, the characteristics {X, V) are close to the characteristics 
{X, V). To that purpose we first prove that (X , V) is close to {X, V), and that {X, V) is close 
to {X, V) (as c +00). 



1. Using the approximation Lemma 14.31 and defining B and E by (|4.29p . we can choose c 
large enough (larger than some ci(T) > 0), so that the characteristics {X,V) are arbitrarily 
close to the characteristics {X,V), and hence satisfy (|4.36p . 
Indeed, using (j4.27p - (|4.28p . we can set 



d:=E-Vlp-E + 



B-B 



Bq dx 



and observe that its norm can be made arbritrarily small, for c large enough. Thus that the 
characteristics {X,V) satisfy ()4.36p follows from the definition of (/, V^) and an application 
of Lemma 14.21 



2. It remains to prove that, the characteristics [X,V) are arbitrarily close to [X^V) if c is 
large enough (larger than some C2{T) > ci(T)). The convergences of E and B to zero as c 
tends to infinity are only weak in time, but, as we will see, this is sufficient to establish this 
property on the characteristics. 
Let us introduce 

{Y,W) := {X,V)-{X,V). 



Then {Y, W) satisfies 
dV 



dt 



dW — — ~ V^ — 
— = E{X) - E{X) + —B{X) 



—BiX) - EiX) - —B(X). 
c c 



We remark that the right hand side of the equation of W is bounded uniformly with respect 
to c. In the same way, we notice that X is Lipschitz with respect to t, uniformly in c and that 
V >—?■ V is Lipschitzian with a Lipschitz constant independent of c. Multiplying the equations 
by Y and Vl^ by y and W respectively and integrating over time, one deduces that for some 
constant C > independent of c, one has 



d 

(4.37) —iY^ + W^)<C{Y^ + W^) + 

dt ,iQ 



-EiX) 



V 



■B{X) -Wds. 



Denote 



V 



Ii := / E(X)-Wds and /a := / B(X)-Wds. 

Jo Jo c 

Let us show that uniformly on [0, T], one has Ii ^ as c +00. Using the definition of E, 
and denoting efc(x) = exp(i27rA; ■ x), we have 



h : -- 



fcez2\{o} 



kAB'^jO) 



sin(sc|- 



+ 



^'(0) - ^'^''(O) cos{sc\k\) 



CkiX) ds. 
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Integrating by parts, we infer 




fc Ai?^(0) 
\k\ 



cos{tc\k\) + [^^'(O) - T^p'^i^)] sin(te|A;|) efc(X). 



cos(sc|A;|) - r^^fO) - ^/(0)1 sin(sc|A:|) ds 



Using dsiWckiX)) = dsWek{X) + W{dsX ■ k)i2Trek and the regularity of B{0), E{0) and 
p{0), we obtain the claim that /i — > as c ^ +cx3. The proof for I2 is similar and therefore 
omited. 

Using the convergences of Ii and I2, ()4.37p and Gronwall's lemma yields the claim on the 
characteristics. Therefore ()4.36p is satisfied for c large enough. 



4.4. With an arbitrary control set? The strip assumption is only used so that the local 
conservation of charge and the zero-mean current property are satisfied by the reference solu- 
tion for Vlasov-Poisson, so that one can apply the approximation Lemma [4.31 More precisely, 
this geometrical hypothesis is only used to get ()4.15p . which in turn allows to ensure these 
properties. For the case of an arbitrary open control set uj, keeping the notations which follow 
Lemma |4.1| we would need to construct a function v satisfying the stronger property: 



where T>' is a line orthogonal to V and which does not meet uj (reduce oj if necessary). 
Unfortunately, we were not able to build such a function. Maybe the use of elliptic functions 
could be helpful in this general case. 



To conclude the proof of Theorems [H [2] and [31 we can now rely on the reference solutions 
which were built in the previous sections to build a solution of the Vlasov-Maxwell system, 
which takes into account the initial data, which remains "close" to the reference solutions and 
finally reaches the state 0. 

As the completions of the proofs are now very similar, we will focus mainly on the construc- 
tion for Theorem [2] (and for Theorem [3] which is a corollary of Theorem [2]). We will ultimately 
explain the needed modifications for Theorem [TJ 

We will propose a functional framework which is bit different to those proposed for the case 
of the Vlasov-Poisson system with external magnetic field in |14j and which is inspired by the 
work of Asano fT] and Wollman [29l[30] . 

This section is organied as follows. We start by proving Theorem [2l This is done in three 
steps. First, we construct a solution of the Vlasov-Maxwell system with source (together with 
the source itself) by a particular fixed point scheme (Subsection 15. ip . Then we prove that 
this solution is relevant (Subsection 15. 2p : this proves Theorem [2] when /i = outside u. The 
general case /i 7^ is explained in Subsection 15.31 Next in Subsection 15 . 41 we deduce Theorem 



□ 




5. Construction of a solution of the non-linear problem 
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[31 In Subsection 15. 5[ we prove Theorem [TJ The last subsections propose several possible 
generalizations of these results. 

5.1. Fixed point scheme. We recall that xq & u and tq > are such that B(xq, 2ro) C w. 
We introduce R > such that 

(5.1) Supp fo, Supp /i C t2 X B{0, R/2). 

We apply Proposition 14.21 we fix / the reference function and T the control time given in this 
proposition. We will assume that c > c(T). 

Let e £ (0, 1). We begin by introducing the operator V whose fixed point will give a solution 
to the nonlinear problem. We first define the domain Ss^b. of V by 





9 G CbiQr) 1 




a. 


Vt G [0,r],Vx G T2, Supp„ g{x,-) C 5(0, i?), 




b. 


b - / L°°([0,T];/f^(T2xR2)) < 




c. 


b - / hyl>°°([0,T];fl'2(T2xR2)) < ^> 




d. 


Vt G [0, T], /^2^jg.2 g{i, X, t;) = xk2 /o(^ 


v)d 


e. 


Vt G [0,T],Vx G T^jSf /jg2 5 c^^^ + divx f^2vgdv -- 


= 0, 


f. 


g\t=o = fo\- 





(5.2) 



Given /o small enough and satisfying ()5.ip . one has f + f & Sg^R (denoting by / the solution 
of the relativistic free transport equation with initial data /o), and consequently S^^r ^ 0. 
From now on, this is systematically supposed to be the case. 

As in |13| 114) . we introduce the following subsets of S{xq, 2ro) x M?: 



(5.3) 

(5.4) 

(5.5) 
(5.6) 



7 



{x,v) G 5(xo,2ro) xR^ / \v\> - and v ■ u{x) < 



(x, v) G S{xq, 2ro) X / > 1 and v ■ v{x) < 



1 

W 
1, 

{x,v) G 5(xo,2ro) x / > 2 and v ■ ^{x) < --\v 



7+ := |(x, v) G S{xo, 2ro) x / ^ . 2,(3,) > gj , 



where i/(x) stands for the unit outward normal to the sphere S{xQ,2rQ) at point x. One can 
readily check that 

dist([S(xo, 2ro) X \ 7^"; 7^") > 0. 



satisfying 



We introduce a C°° Pi smooth function U : S{xQ,2rQ) xR"^ ^ 

( 0<U <1, 

(5.7) <^ i7 = 1 in [Sixo, 2ro) x R^] \ 72-, 

[ J7 = in 73-. 

We also introduce a cut-off function T : — ?• IR+, of class C°°, such that 

T 23T 



(5. 



T = in 



T' 




\^7T 1 




U 




"■48. 




[ 48 ' J 



and T = 1 in 



24 24 



26 



OLIVIER GLASS AND DANIEL HAN-KWAN 



Now, given g € Se^R, we associate E^,B^ in 
solution of the Maxwell equations: 



\T];H^{T^))nC\[0,T];H^T^)) as the 



(5.9) 



dtEf 



gvidv, dtEl + cd:,,B3 



gv2 dv, 



dxiEf + 8x2^2= / gdv- gdvdx, 



Ea. B3 



t=0 



Ba. 



E'\t=o 
Then, we define 

(5.10) V{g) := f 

to be the solution of the following linear system 



f f{0,x,v) = fo onT^ X 



(5.11) 



in [0,r] X [(T^ X M^) \7-], 

c 

[ f{t+, X, v) = [l- T{t)]f{t-,x, v) + T{t)U{x, v)f{t-,x, v) on [0, T] x 7" 



dtf + v.V^f + {E3 + —B3).VJ 

c 



In order to explain the last equation of (j5.1ip . we introduce the characteristics {X,V) 
associated with the force field E^ + ^B^. In this equation, f{t~,x,v) is the limit value of / 
on the characteristic {X,V){s,0, x,v) as time s tends to t~ . (Observe that for times smaller 
than t, but close to t, the corresponding characteristic is not in 7~.) When the characteristics 
{X, V) meet 7" at time t, then the value of / at time is fixed according to the last equation 
in ()5.1ip . Loosely speaking, one can interpret the function [1 — T(i)] + T{t)U{x,v) as an 
"opacity" factor which varies according to time, the modulus of the velocity and to the angle 
of incidence of the trajectory on S{xQ,2r()). In this process, a part of / is absorbed on 7", 
which varies from the totality of / to no absorption. Typically, when the velocity and the 
normal incidence of the trajectory are sufficient, the particles are absorbed. We now consider 
a continuous linear extension operator 

7f : H^{[T'^ \ B{xo, 2ro)] x R^;R) H^{T'^ x R^;R), 

which continuously maps H%[T'^ \B{xo,2ro)] x R'^;R) to H^J"^ x R'^;R), for any s G [0,3] 
(by classical results, it is possible to consider such an extension operator because B^xq, 2ro) is 
smooth). In addition, it is clear that if a function / has compact support in v, we can ensure 
that vf(/) also satisfies this property. 

Then we modify this operator in order to get a new operator vr which has the further 
property that for any integrable / S H'^{{T'^ \ B{xq, 2ro)) x M^), one has 



(5.12) 



7f(/) dv dx 



T2xR2 



fo{x, v) dv dx. 



This condition can easily be obtained by considering a regular, compactly supported, nonneg- 
ative function fi with integral 1 over B{xQ,2rQ) x and such that J vfidv = 0, and adding 
to vf(/) the function 



7t2xR2 Jt2 
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We further modify the operator so that it satisfies the local conservation of charge. To this 
end, we add another regular distribution function h, compactly supported in B{xQ,2rQ) x M^, 
such that 

hdv = and div / hv dv = —dt / 7r(/) d-y — div / TT{f)vdv. 



This can be done as in the proof of Proposition 14.11 but here this is much easier to build such 
a function since here we do not need the zero mean current condition (j4.6p . and this can be 
achieved without making any assumption on the control set. We finally obtain a continuous 
affine extension operator vr. 

It is convenient to introduce another time-dependent function T depending on time such 
that: 

T 



(5.13) T = in 



"•100 



and T = 1 in 



— T 

48' 



Next, we introduce the continuous operator 11 given, for / : ([0, T] x [T^ \ B(xo, 2ro)] x M?) U 
([0,r/48] x T2 X M2) by: 

(5.14) (n/)(t, X, v) := (1 - fit)) fit, X, v) + fit)[7Tfit, ; ■)]ix, v). 
We finally define V[g] for g G S^^r by: 

(5.15) V[g] := f + '^if\{[o^T]x[T2\B{xo,2ro)]xR^)U{[0,T/i8]xT2xR2)) m [0, T] X X M^. 

The goal is now to prove the existence of a fixed point for this operator. For clarity of 
exposure, we cut the argument into several lemmas. 

Lemma 5.1. The set Ss^r, endowed with the topology, is a convex compact subset of 
CfUQx), and V is continuous in the topology. 

Proof of Lemma \5. 1\ That is a convex subset of Cj^iQ-r) is straightforward. That it is 
compact can be proved exactly as in |13| Section 3.3] or |14) . relying on Ascoli's theorem. 

Likewise, the continuity can be proved as in |13| Section 3.3] or |14) and we therefore omit 
the proof. Note that by Sobolev imbeddings, and are Lipschitz-continuous in x. □ 

Lemma 5.2. There exist Eq > and Rq > such that for any < e < Eq and R > Rq, we 
have: 

ViSe,R) C S,,R. 

Proof of Lemma\5^ Let g € S^^r'- let us prove that Vig) G Sefi- That d., e. and f. hold is 
true by construction of V. 

That a. and b. hold comes from the fact that one can first prove: 

a'. Vt G [0,T],Vx G T2 \S(xo,2ro), Supp^ V(5)(t, x, •) C 5(0, i?), 

b'- l|1^(5)llLoo([o^r];H3(T2\B(a:o,2ro)xIR2)) < 

as soon as e is small enough and R large enough. 

These estimates are obtained with the same method as in the classical proof by WoUman 
|30) of the existence of solutions of the Vlasov-Maxwell system (without any source); special 
emphasis will be put on the new terms coming from the absorption procedure on Six^, 2ro), 
that actually cause significant technical difficulties. The principle is to perform energy 
estimates separately on [T^ \ ^(xo, 2ro)] x and on BixQ, 2ro) x (instead of on x 
as in Wollman's paper), by using the Vlasov equation satisfied by Vig), and then make the 
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sum of the two contributions. The reason for this treatment is the discontinuity across the 
hypersurface S'(xo,2ro) which is caused by the absorption procedure. 

In order to simphfy the notations, we denote by hi the function V{g) restricted to [T^ \ 
B{xo,2ro)] X and by /i2 the function V{g) restricted to B{xo,2ro) x M^. 

The estimates, first in L^, then in H^, etc., are obtained by induction. Let us start with 
the L norm. By multiplying by hi the equation satisfied by hi and integrating over 
-B(rEo,2ro)] x M?, we obtain after an integration by parts: 



d 



\hi\'^dxdv= I \hi'^v ■ i>{x) da dv. 

;T2\B(xo,2ro)]xR2 Js(xo,2ro)xR2 



Similarly, we obtain for /i2: 



d_ 
IE 



I \h2'^dxdv = — / \h2\^v ■ v{x)da dv. 



For clarity, it is convenient to introduce the following notations: 

7™ := {{x,v) e S(xo,2ro) xR^, v {-u{x)) > 0}, 
7°"* := {(x, v) € S{xo, 2ro) x M^, v (-z^(x)) < 0}. 

Using the definition of the absorption (see the last equation of (|5.1ip ). we have 
(5.16) / \hi\'^v ■ {—^{x)) da dv + / Ihi]"^!) ■ { — ^{x)) da dv 

J -yin J yout 

\h2\'^v ■ {—u{x)) da dv — j \h2\'^v ■ {—u{x)) da dv 

,out 

X [1 — (i, X, v)]v ■ {—v{x)) da dv, 

for some smooth function ai satisfying < ai{t,x,v) < 1. Therefore the right hand side of 
(|5.16p is non-positive and one obtains estimates for V{g): 



d_ 
di 



\hi\'^dxdv+ / |/i2pdx(if < 0. 

T2\B{a;o,2ro)]xR2 J B{xo,2ro)xR'^ J 



Thus, using (jS.lSp . one obtains estimates for V{g): 

[ \Vig)\^dxdv<C\\fo\\mT2y,^2), Vte[0,T]. 



T2xR2 



Let us now treat the norm. We first differentiate the Vlasov equation with respect to 
space or velocity and we obtain as before (here D := dx^ or d^-, i = 1, 2): 



(5.17) ^[ _ \Dhi\'^dxdv= [ \Dhi\'^v u{x) da dv 

dt J[T2\B(xo,2ro)]xR2 JS{a;o,2ro)xR2 

- / _ Dv-V^hiDhidxdv- [ _ D{E3 + —3^) -V^hiDhi dx dv. 

ifT2\B(xo,2rn)lxR2 y[T2\B(a;n,2rn)l xR2 C 
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Similarly, for /12, we obtain the following estimate: 

(5.18) ^[ \Dh2\^ dx dv = - f \Dh2\^v ■ v{x) da dv 

- [ Dv- V^h2Dh2 dxdv- [ D{E^ + —B^) ■ VMDh2 dx dv. 

One can bound the last two terms of equations (|5.17p and ()5.18p (using the Cauchy-Schwarz 
inequality) by: 

C(l + ||(i?^i?^)|Uj.(^3))(||/^l||2,^ + ||/i2|||/J. 

It remains to estimate the first terms in the right hand sides of (|5.17p and (jS.lSp . One can use 
the definition of the absorption procedure in ()5.1ip to compute the derivatives dth2, dyh2 and 
57-/12 on 7*" (denoting by r a tangential unit vector on S{xq, 2ro)) in terms of those of hi on 
7*"'. To obtain the "missing" derivative 5y/i2, one uses the Vlasov equation, for {x^v) € 7*". 
We can obtain the following relation on 7*": 

(5.19) Vt,x,vh2{t,x,v) = a{t,x,v)Vt,x,vhi{t,x,v) + l3{t,x,v,E^ ,B^)hi{t,x,v), 

for t G [0,r],x G S{xo,2rQ),v G M^ 

for some smooth bounded functions a and /3 satisfying < a < 1. We refer to Section 
3.3] for full details. 

Summing the two contributions of Dhi and Dh2, we deduce that the boundary terms can 
be rewritten using ()5.19p : 



(5.20) / \Dhi\^v-{-v{x))dadv+ I \Dhi\^v ■ {-v{x)) da dv 

\Dh2\^v • {—iy{x)) da dv — / \Dh2\'^v • (— z^(x)) da dv 

J ^out 

\Dhi\^ X [1 - a^{t,x,v)]v ■ {-v{x)) dadv 

P'^\hi'\^v ■ {—v[x)) dv — i aPDhihiv ■ [ — v{x)) da dv. 



The first term in the right hand side is non-negative. As we will see, the second one will 
be absorbed using a Gronwall estimate. We have to handle the last term, which is the most 
difficult one, in order to write it with the same form as the second one. We first claim that 
on 7*", we can write (with again D := dx^ or 9^,. , i = 1, 2): 

Dhi = ai{t, X, v)dthi + 02(^5 x, v)drhi + a3(t, x, v)dvhi, 

where the are smooth functions (depending of course on the partial derivative D). This is 
simply proved by using the fact that hi is the solution of a Vlasov equation. Then we can 
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compute the following integral: 

t r rt r /j2 

aPDhihii) ■ {—u{x)) da dv ds = / / al3ai{t,x,v)dt—v ■ {—i>{x)) da dv ds 



JO J7 



2 

2 



+ / a(3a2{t,x,v)dr—v ■ {—iy{x)) da dv ds 
Jo J-r^" 2 



/j2 

+ / / af5a-i{t^x^v)dy-^v ■ {—i>{x)) da dv ds. 



'0 J 

At this point, we can integrate by parts with respect to i, r and which yields only two 
boundary terms in time (at time and t respectively, coming from the integration by parts in 
time) and terms of the form: 

/ / (i2{SiX,v)\hi'\^v ■ { — v{x)) da dv ds^ 
Jo Jy'" 

for some smooth bounded function ^2- The boundary terms in time are the two following one: 



ai 



=ol3\s=ofoV ■ {-iy{x))dadv, 



and 



r 



\s=tf3\s=thi\s=tv ■ {-'^(x)) da dv. 



To treat the second term, we use a trace estimate and Young's inequality, as follows: 
/ af3h\v ■ {—v{x)) dxdv < C\\hi\\^2\\hi\\jji 

<cm\i. + ^\wh,\\i,. 

Using this computation we obtain: 

J V I / _ \Dhi{t)\^ dxdv+ I \Dh2{t)\^ dx dv\ 

<C\\h\\L^^C V I I _ \Dhxf dxdvds + CWhxf^^ 

r.^<aa l ^0 J [T2\B(xn,2rn)l XR2 



+ V ( / _ \Dfo\^dxdv+ [ \Dfo\^dxdv], 

and we conclude by Gronwall's lemma. 

The higher order derivatives are obtained in the same way by induction. In particular, we 
observe that we can obtain a similar formula as ()5.19p for higher derivatives. Indeed, let us 
prove that for i,j G {t, a^i, a;2, ^2) "^2} we have: 

(5.21) dljh2 = adijhi + ^Ijdjhi + ^Ijd.hi + 7^ j/ii on [0, T] x S{xo, 2ro) x 

where a is the same function as in ()5.19p and f3jj,f3fj,'yij are smooth bounded functions of 
(t, X, v). From ()5.19p . we obtain easily that ()5.2ip is true for any j when i is a tangent direction 
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to [0,T] X S{xQ,2rQ) x M?. For the remaining (normal) direction, we differentiate the Vlasov 
equation satisfied by hi and /i2: 

dtidjhk) + V ■ V^idjhk) + F3 ■ V^idjhk) = - [djv ■ V^hk + djF^ ■ V^hk] , 

denoting by := + —B^ and k = 1,2. Then we replace the first order derivatives of /12 
using (|5.19p and the second order derivatives using (|5.2ip except the one involving the normal 
direction in v ■ Va;(c?j/i2), that is [v ■ V x{djh2)- As a result, we can express this last term 
in terms oi a{v ■ v)v ■ Vxidjhi) and lower order terms as in ()5.2ip . 

Then, when one studies the contributions from 7" in the estimates, the most harmful 
terms come as before from the crossed terms adijhi/Sijdjhi, but as before, one can use an 
integration by parts argument. Similar statements hold for derivatives of order 3. Note that 
when proving the and the estimates, the time derivatives have to be understood in the 
weak sense. 

Once a', and b'. are obtained, that a., b. hold for V{g) is a consequence of the construction 
of the extension operator. 

There remains to prove that c. also holds. This is actually a consequence of b'. and the 
fact that hi solves the Vlasov equation, so that 

Pt^l|lLoc([o,T];H2([T2\B{xo,2ro)]xR2)) - 

Finally that c. holds is obtained by construction of the extension operator. 

□ 

We conclude this section with the following statement. 
Lemma 5.3. If R > Rq and < e < eo, V has a fixed point in Ss^r. 

Proof of Lemma \5.S[ This is an immediate consequence of Lemma |5.H Lemma 15.21 and the 
Schauder fixed point theorem. □ 

5.2. Relevancy of the solution. Let us now prove that the solution we have just built is 
relevant for the controllability problem of Theorem [21 

Lemma 5.4. There exists ei > such that for any < e < ei, all the characteristics (X, V) 
associated with + , where f is a fixed point ofV in Ss^r, meet for some time in 

U2' 12 J- 

Proof of Lemma \5.4\ This is a consequence of Proposition 14.21 and the fact that a fixed point 
(/, E^ , B^) is a small perturbation of the reference solution (/, B). When e is chosen small 
enough, we have: 

||S^-^||oo<l, 

\\Ef -VTpW^^l. 

This involves that the characteristics {X, V) associated with the force field E^ + —B^ still 

satisfy the property: 

(5.22) 

Vx G T^.Vv G M^3^ G [r/ll,10T/ll], X{t,x,v) G S(xo,3ro/2), and \V\{t,x,v) > 3. 

As in |14) . we can prove that this implies that, provided that e@ is small enough, the trajec- 
tories meet 7 during [j^, ^^], so that all particles are absorbed. The technical points are 
almost identical to those of the proof given in Lemma 5.3], and so we omit them. □ 
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We can finally end the proof of Theorem [2] when /i = outside uj. 

Proof of Theorem when /i = outside u . By Lemma 15.31 and Lemma 15. 4^ we obtain a so- 
lution / such that: 

Supp[/(T,.,-)] CwxM2. 

We obtain furthermore that the distribution function / satisfies for some function G supported 
in uj: 



dtf + v ■ V^/ + div.„ 



(Ef + -v^Bf)f 
c 



G in [0, T] X X 



The function G is then our control function. Steering finally the state exactly to /i can be 
performed via a function supported inside tj, which also yields a source term internal to uj. 
This ends the proof of Theorem [2] (in the case where the target /i = outside oo). □ 

5.3. From /i = to general targets. Up to now, the proof of Theorem [2] has been com- 
pleted for the case /i = outside uj. We cannot directly invoke a reversibility argument to 
get the general case, since we cannot precisely say at which state the electro-magnetic field 
[E, B) is at final time. 

Nevertheless, it is possible to modify the fixed point scheme to complete the proof in the 
general case. Keeping the same notations, we introduce a function / which is the solution of 
the Vlasov equation: 

' f{Q,x,v)=fi onT^xM?, 
(5.23) I aJ + 7)-Vj+(i^9 + i^ij9).v,/ = 0in [0,r] X [(T2 xR2)\^-]^ 

^ f{t+,x,v) = [I - T{t)]f{r ,x,v) + T{t)U{x,v)f{r ,x,v) on [0,T] X7-. 
We observe, as in the proof of Theorem [2] (using the property of the characteristics), that 
f\t=T = 0. 

We finally define the fixed point operator V which maps g ^ r to: 
(5.24) 

Vb] := / + n((/ + /(r-t,X,-u))|([o,T]x[T2\B(xo,2ro)]xR2)u([0,T/48]xT2xI82)) ^ [0, T] X X M^, 
and conclude as previously. 

5.4. Proof of Theorem [3l As already explained. Theorem [3] is a consequence of Theorem [2j 
We recall that the speed of light here is equal to 1 . Let us show how to reduce to a sufficiently 
large c. The principle is to rely on the scale invariance of the equations. Indeed, taking A 
large enough in (|2.5p increases the value of the "effective" speed of light. 

Let us be more specific. Given a; (satisfying the strip assumption) and 60, we apply Theorem 
[21 We obtain a time Tq, and apply the result to T = Tq. We obtain a constant ac > and a 
speed of light ctq > 0, such that whenever c > ctq, starting from /o,-Eo,/i satisfying (jl.lUp - 
()1.13p . we have a solution of Vlasov-Maxwell with speed of light c, steering /o to /i in time 
To. 

We define Ai := I/ctq. For A G (0, Ai], we observe that 

c := 1/A 

satisfies c > ctq. We define ki := min(l,A^)K;. Given /o,/i,ii^o satisfying (|1.14p - (jl.l5p and 
such that 

II(/o,-^o)||h3 < II/iIIh3 < 



ON THE CONTROLLABILITY OF THE RELATIVISTIC VLASOV-MAXWELL SYSTEM 



33 



then 

fo{x,v) := fo{x,Xv), fi{x,v) := fi{x,Xv), Eo{x) := l/A^£^o(x) 

also satisfy (|1.10p - ()l.lip and moreover (I1.12p . We also introduce the initial magnetic field Bq 
(and its rescaled counterpart Bq): 

Bo{x) := Xboix), Bo{x) := l/X^Bo{x). 

Therefore, we observe that: 

Bo{x) = cbo{x) 

and we can apply Theorem [2] to obtain a control function Gx such that the solution of the 
system: 



dtf" + vx ■ V,/-^ + div, (^^ + XviB^)f^ 



dtE^ + -curlS^ 
A 



f^vxdv, 



dtB^ + -cm\E^ 
X 



div E 



A 



/|t=0 



with Vx 



(5.25) 



satisfies: 



f^dv 

2 

^|t=0 = 



rdvdx, divB^ = 0, 



En 



Bq (= cbo 



1+ 



f\t=To ~ /l- 

One deduces a solution (/, i?) to the Vlasov-Maxwell in the original variables by performing 
the "inverse" change of variables: 

f{t/X, X, Xv) := f^{t, X, v), E{t/X, x) := X^E^{t, x), B{t/X, x) := X^B^it, x). 

In particular, we observe that 

f\t=To/X = /l- 

which gives the controllability in large time of Theorem |3l 

5.5. Proof of Theorem [Tl In order to prove Theorem [H we modify the previous construction 
as follows. First, we recall that in the construction of the reference solution in Section [3l we 
have considered the sequence of sets uj' C oo" C oo" C OJ. Let / be the reference solution given 
by Proposition 13.11 

Let us describe the principle of the absorption in this case. We consider a covering of balls 
{B(x, rQ))x£uj" of oj", whose radius tq is small enough so that the covering with balls of radius 
2ro is included in uj. By compactness, we can assume that this covering is made of a finite 
number of balls B{xi, vq), i?(xfc, tq). 

On each ball B{xi, 2ro) (for i = 1, ...,k) we impose the same absorption procedure as before. 
Then one can proceed with almost the same proof as Theorem [2] The only difference appears 
in the proof of Lemma [5.21 we now have to estimate the norms L^, H^, etc, in each connected 
component of \ U^^^5(xfc, 2ro). Remark that there is only a finite number of such regions. 
Then, we make the sum of all contributions of each region to get the relevant estimates. 

Finally, we obtain the construction of a control supported in u which brings the state / of 
the system to a value satisfying: 

Supp[/(r,.,-)] Cu;xM2, 
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for T large enough. Taking T > T large enough so that one can apply Theorem H] (control of 
linear Maxwell equations), we can further modify this solution so that 

Supp[/(t, •, •)] C w X for any t G [T, f], 

and so that at time T, we have 

/(T,-,.)=0 on{T^\u)xM.\ 

and 

{Eif,-),B{f,-)) = {0,0). 

By a slight abuse of notation, we still write T instead of T. We clearly obtain that the 
distribution function / satisfies for some function G supported in cj: 



dtf + v.Vxf + div^ 



(^/ + -v^Bf)f 
c 



G in [0, T] X [T^ \uj]x 



The function G is then our control function, which brings the state to (0,0,0). 

We can finally directly invoke a reversibility argument (take A = — 1 in (|2.8p ) or proceed as 
in Subsection 15.31 to steer the system to general (fi, Ei, Bi). 



5.6. A C or Cauchy theory ? In the Vlasov-Poisson case, in order to build the solution 
in a framework, we used the same absorption procedure which roughly consists in absorbing 
fast particles which enter the control zone with sufficient normal incidence. This procedure 
is a bit costly in terms of regularity, and it causes a loss of moments in velocity. In \13\ I14j , 
this is compensated via interpolation by a loss of regularity. Nevertheless, this regularity is 
re-gained afterwards using the elliptic regularity of the Poisson equation. One first difficulty 
in the Maxwell case lies in the fact that this is not possible to do so (as the Maxwell equations 
are hyperbolic) . A more severe obstruction lies in the hypothetical application of the Glassey- 
Strauss-Schaeffer Cauchy theory for Vlasov-Maxwell. Since the Maxwell equation can be 
seen as a wave equation, it is quite remarkable to be able to build solutions in a framework 
(even worse, in the wave equation form, there seems to be a loss of derivative). Actually, one 
can observe that the theorems of Glassey-Strauss-Schaeffer crucially rely on the fact that the 
distribution function is a solution to a Vlasov equation (without source). This allows to trade 
derivatives like {dt + v ■ V^) for derivatives in v. Then, integrations by parts in v allow these 
authors to perform suitable estimates. For our problem, this is unfortunately impossible 
to do so, since in the control zone, the Vlasov will feature a source (the control function) with 
only regularity. Therefore it seems difficult to build solutions with these methods. 

Instead of the theory, we have thus opted for theory of Asano and Wollman (although 
it is less sharp in terms of regularity). 



5.7. Removing the smallness assumption on B in Theorem [H Let us explain how one 
can remove the smallness assumption on Bq (and Bi) in Theorem [TJ When one applies the 
return method, instead of building a reference solution which goes from (0,0,0) to (0,0,0), 
one builds instead a solution going from (0,0, -Bj) to (0,0,0). To do so, one can still rely on 
the controllability of the Maxwell equations. Then we apply the same procedure as above. 
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5.8. Generalization to data with unbounded support in velocity in Theorem [H 

Until now, we have dealt with distribution functions having a bounded support in velocity. As 
claimed before, we expect the theorems contained in this paper to be generalized to distribution 
functions with unbounded support in velocity, but having exponential decay at infinity. 

To that purpose, instead of the standard Sobolev space of the local theory of Wollman 
|ljU) . we may use weighted Sobolev spaces (accounting for the exponential decay in velocity), 
following the work of Asano [1]. It is likely (computations get more tedious) that the fixed 
point procedure proposed in this section is still relevant. 

6. Appendix: Proof of Lemma 14.21 

In this section, we provide a proof of Lemma 14.21 Actually we establish a result which 
is more general, by allowing more generic magnetic fields. Let us first describe the class of 
magnetic fields b that are admissible. 



Definition 6.1. // a magnetic field b is such that either b or —b satisfies the conditions 1. 
and 2. below, we say that b satisfies the bending condition. 

Condition 1. Geometric control condition on b. We assume that there exists a compact 
set K of on which b > and which satisfies the geometric control condition: 

(6.1) Vx e Ve e §^ 3ye M+ such that x + yeeK. 
For a compact subset K of and r > we denote 

(6.2) Kr := {x G TV d{x, K) < r}. 

The first geometric assumption can be reinterpreted with the help of the following lemma (see 
|14) for an elementary proof). 



Lemma 6.1. Let K cT"^ such that b > on K and satisfying ()6.ip . Then there exists b>0, 
d > and D > such that 



(6.3) b>b onK2d, 

d' 

t,t + - 
' 2 



(6.4) Vx G Ve G §\ 3t G [0,L>], Vs € 



X + se G Kfi- 



Roughly speaking, this lemma gives, for a ray of light (whose velocity is a unit vector), the 
maximum time it can spend outside during one passage (that corresponds to D) and the 
minimal time it has to spend in K^, during one passage (that corresponds to d/2). 

This allows us to introduce the second condition on the magnetic field: 
Condition 2. Bound from below. We assume that there exists b G M" such that 

b > 6 on T^ 

and such that, keeping the same notations as in Lemma |6.H 

(6.5) bD + b^ >0. 

It is clear that a constant nonzero magnetic field satisfies the bending condition. 

A less general notion was introduced in |14) (the second condition was replaced by a fixed 
sign assumption, which corresponds to 6 = 0); the main interest of this new definition comes 
from the fact that it is stable by small perturbations. Indeed, if some b satisfies the bending 
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condition, then there exists Eq > such that any b with \\b — b||oo < ^O) stih satisfies the 
bending condition. In particular, the following lemma allows to improve some of the results 

of [n]. 

Lemma 6.2. Let xq € T'^,rQ > and M > 1. Given b € C^(T^) satisfying the bending 
condition, there exist: 

• cq > depending on b,Xo,ro, 

• m > depending only on b,Xo,rQ, 

• T > depending on b, xq, and M , and 

• K depending on b, Xo,ro and M, 

such that for all ^ e L°°(0, T; ^^^'"^(T^ x M?)) satisfying ||5||l°° < k, if c > cq then the 
characteristics {X,V) associated with i)^b +5 satisfy: 

(6.6) Vx e T^Vt; G such thatM> \v\ > m, 3t G {T/4,3T/4), X{t,0,x,v) G B{xo,ro/2) 

and for all s G [0,r], y < \V{s,0,x,v)\ < 2\v\. 

We prove Lemma [6. 21 in several cases of increasing complexity. In a first time (Cases 1-3), 
we suppose that 5 = 0. In Case 4, we explain how to take the additional small force 5 into 
account. 

In all cases, we define 

(6.7) b := max b(x). 
Before giving the proof, we give its main ideas. 

Idea of the proof of Lemma 1 6. ^ The proof is a variant of that given for [141 Proposition 
5.1]; here we have to be careful of the effects of relativistic transport and in addition, the 
assumptions on b are more general. Conversely, the assumption of smallness on 5 will allow 
us to simplify some aspects of the proof. 

In a first time, we consider non-relativistic trajectories (c = +oo) and ^ = 0. One can 
show that for (non-relativistic) free transport (whose trajectories are straight lines), there is 
only a finite number of bad initial directions which produce trajectories which never meet the 
ball B{xQ,rQ). Indeed, trajectories with directions having irrational slopes are dense in the 
torus and therefore have to meet B(xo,ro). Considering periodic trajectories with rational 
slopes, because of the "thickness" of the ball B{xo,rQ), only a finite number of them produce 
trajectories which do not meet the ball. A magnetic field satisfying the bending condition 
allows to circumvent these bad initial directions. The proof of Lemma 16.21 is based on the 
three following ingredients. 

• Let us define the classical characteristics {X,V) associated with v^b. We can first 
observe that when the modulus of the initial velocity \v\ is large, the characteristics 
(X, V) remain close to those of free transport (straight lines) at least for a time of 
order l/\v\. This allows to show that trajectories X which initially have a velocity 
with a good direction, meet the control zone as it is the case for free transport. 

• For trajectories which initially have a bad direction, we rely on the "rotation" effect 
provided by the magnetic field, which yields that bad directions become good directions 
after some time r. The condition 1. ensures that trajectories "often" meet the zone 
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where b > 0, in which the direction of the velocity gets rotated and therefore do not 
remain of bad direction. Note that at each passage in the zone b > 0, the angle of 
the velocity is only modified by ©(l/lul), but using the reinterpretation of Lemma 
16.11 we can see that during a time of order 1, there are at least 0(1^1) passages in 
this zone. On the other hand, the condition 2. in the bending condition allows to 
say the trajectories are not too affected by their passages in the zones where b can be 
non-positive (and which could otherwise annihilate any bending effect). 
• Finally, by choosing c sufficiently large and 5 sufficiently small, by a perturbation 
argument, we get that the relativistic trajectories {X, V) are close to the classical 
trajectories (X,V), which allows to conclude. 




{b > 0} 

^ 



Figure 2. An illustration of the proof of Lemma 16.21 the bending effect of 
the magnetic field. 



Proof of Lemma\6i 

In the sequel, we denote (X*,y*) the characteristics associated with the relativistic free 
transport and (X, V) the characteristics associated with relativistic transport with the force 
v^h + ^. 

Case 1. Constant magnetic field and J = 0. Let us first suppose that b is constant and 5^ = 0; 
for readability we assume here that b := 1. 

As noticed in |13| Appendix A, p. 373-374], there is only a finite number of directions in 
S"*^ for which there exists a half-line in which does not intersect i?(a;o, ^q/S). Identifying 
with [0, 27r[, we denote them ai, ...,ajv G [0, 27r[. 

We introduce the neighborhoods of Oj: 

Vi = (a, -/3/2,Q, + /3/2), 
as follows. Let /3 > and r > satisfying: 



(6.8) /3 < mind(Vi,V,)/8 and r = . / 1 + ^ min d(Vi, V,)/7. 

By construction, observe that: 

T 

/5< 



^0 

r2 



By a compactness argument, there exists L > such that for any x £ T , and for any 
Oj € \ U^^^Vj, any trajectory X* starting from x with a direction has to cover at most 
a distance L to meet i?(xo,ro/8). 
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We choose cq and m large enough such that: 



, , La/1 + m2/c2 / 1 1 

(6-9) := ^ ^ < Z3 + 3 < ^- 



This is the time "free" trajectories X* with velocity of modulus m take to cover the distance 
L. We observe that for any \v\ > m, we have 



(6.10) Tui := < Tm- 

Now let X € T^jV G with \v\ > m. First note that for all t, \V{t,0,x,v)\ = \v\. Now let us 
discuss according to the direction of v. 

• First case: ^ G \ U^^Vi. 
We have, for any t < T^^^ < Tm, 

\v\ ^'i 1 L\l + \v\yc^) 



\X*{t,0,x,v) -X{t,0,x,v)\ < 



y^T+W/^ 2 ^1 + |t;|2/c2 2\v\ 
2 



< — ^/l/c2 + l/m2. 



We can impose cq and m large enough such that: 

(6.11) — ^l/c2 + l/m2 <ro/8. 

As a result, we obtain that 

3te (0,r„,],X(t,0,x,t;) e B(xo,ro/4), 
and (|6.6p is trivial here since |l/(t, 0, x, f )| is conserved. 

• Second case: ^ G Vjf^^i, say Vj. 

As already explained in the heuristics, we wait for a time r, that is, we consider: 

{x\v') := (X(t,0,2;,7;),F(t,0,x,7;)). 

We observe that because of the "rotation" induced by the magnetic field, the angle between v 
and v' is equal to r/y^l + |f |2/c2, which is larger than t/{\J1 + M'^/cq). 

Consequently, due to the choice of /3, 

and thus we are in the same case as before. 

Therefore, defining T := 4:{Tm + t), we have proven that (actually up to a harmless trans- 
lation in time): 

3t e (r/4,3r/4], X{t,0,x,v) e S(xo,ro/4), 
together with the condition ()6.6p on the velocity field. 
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Case 2. Positive magnetic field modulus and ^ = 0. Here we suppose that b = b(x) > on 
We are in the case where in Lemma [6.11 we can take K = Kd = and 



6 = inf b > 0. 



Keeping the same notations as before, we set r G (0, T] and /3 > in order that 



'1 + ^ 

(6.12) /3 < min(i(Vi,Vj)/8 and t = ^ — min (i(Vi, Vj)/7. 

The proof is very similar to the previous one. Indeed, the fohowing comparison estimate stih 
holds for t < Tj„| < r where T]^| was defined in (|6.10p : 

(6.13) \X*{t,0,x,v)- X{t,0,x,v)\ <6— ^l/c2 + l/m2. 

Let X S T'^,v € M^. We distinguish as before between two possibilities. Using the inequality 
()6.13p . the first case holds identically for m large. For the second case, we just have to check 
that with this magnetic field, the velocity is rotated by an angle at least equal to /3 after some 
time t < T. 

We use the following computation for general {x,v). Denote by 9{t) the angle between 
and V{t,0, x,v). We compute the scalar product of V{t,0,x,v) and £[]£(^^£i£!l ^ using the 
identity: 

dV{t,0,x,v) V{t,{),x,v)^ 
— b(A (t, 0, X, v)) 



'^^ '1 + |F(t,0,x,?;)|Vc2 



We straightforwardly obtain that |y(f, 0, x, f)! = \v\. Then, taking the scalar product of 
and — ^ , we likewise obtain: 

y 1 + \v\'^ /& 

(even if sm6{t) = in which case one considers the scalar product with v.) We deduce that 

e'{t) > h/^i + NP/ci 

Thus going back to (^,'5), by the intermediate value theorem and the definition of the 
neighborhoods Vj and to ()6.12p . there is a positive time t less or equal to r for which we have: 

F(t,0,x,^) GSi\uiIiVi, 

and we conclude as in Case 1. 



Case 3. Magnetic field satisfying the bending condition and 5 = 0. Let us now consider the 
general case for , but still without the additional force ^. 

Given K satisfying the geometric condition ()6.ip . we introduce d and D as in Lemma [6.11 
We consider as in ()6.2p . Using the condition 2. (bound from below), we have: 

(6.15) ^ — Db+H >0. 
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(we recall that by definition, 6 < 0). We assume here that in addition to 
/3 > are such that 

r7 



r € (0, T] and 



2y/l + KP/ci{D + d/2) 



Let X G and e M^. Once again we distinguish between the two possibilities (good or 
bad direction). As before the first case is still similar since ()6.13p is still valid. We have to 
give a new argument only for the second case. 

To that purpose, we examine the behavior of the characteristics during the time interval 
[0,r]. 

• By (|6.4|) . each passage in \ of X*(t, 0, x, v) lasts at most during a time equal to: 



D^/l + \v\'^/c^/\v\. 

The characteristics X are not straight lines since they are modified by the magnetic field. 
Let us prove nevertheless that if \v\ and c are large enough, then the particle X(t, 0,x,u) 
can remain at most during a time D ^J\ + \v^^^/(P' /\v\ in \ Kd- Let x G \ Kd, and 
1^ € S"^, let o" > 0. By Lemma 16.11 there exists some s < Dy^l + |up/c^/|t;| such that 
X#((T + s, a, X, v) € K. Now we can evaluate as before: 



X*{a + s,a,x,v) - X{a + s,a,x,v) < h^yjl/cl + l/|up. 

Hence, we can impose cq and m large enough such that when \v\ > m, 

X{a + s,a,x,v) G Kd- 

Hence, we deduce that each passage of X{t,0, x,v) in \ Kd lasts at most during a time 
equal to Dy^l + |t;p/c^/|u|, which proves the claim. 

• Likewise, given x G \ Kd, £ and cr > 0, by Lemma 16. H there exists some s < 
Dy/l + \v\^/c'^/\v\ such that X#(r, cj, x, G K for all r G [s + a, s + a + ^^Jl + \v\'^/c /\v\]. 
We have for ah r G [s + (t,s + a + ^^Jl + |wp/c^]/|u|] that 

1^1 ((i^ + f)Vi + IHVcVH))' 



|X*(r, fj, X, t;) — X{t, o", x, < h- 



Hence we can choose m and cq large enough such that for any \v\ > m and when c > cq, 
X(r, cr, X, G Kd for ah r G [s + cr, s + cr + f y^l + \v\'^/c\. 

• We see that a trajectory alternates between passages in Kd and in \ Kd- We denote as 
previously by 6{t) the angle between y(t,0, x,f) and v and aim at finding a lower bound for 
9{t). We introduce 



N :-- 



T\V\ 



,/l+W/^{D + d/2) 
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(here [-J denotes the floor function) and 

A := |t € [0,r] j X(t,0,x,t;) G T^\K^. 

There are at most distinct intervals of length ^1 + |up/c^(L' + (i/2) in [0,r]. Using the 
considerations above, one sees that the Lebesgue measure of A satisfies 



\A\ <{N + l)D^/l + \v\^^/c^^/\v\. 
Indeed, taking to := inf A, one has 



\An[to,to + vT+W/^iD + d/2)]\ < D^/l + \v\ycy\v\, 



and we can make the same reasoning with ti := inf(^ H {t > to + y^l + \v\'^/c'^{D + d/2)}, 
etc. Likewise, 

d 
2 

and 



\{[0,t]\A) n[to,to + ^/TTWJ^iD + d/2)]\ > -^l + \v\ycy\v\ 



.d 



\[Q,t]\A\>N-^/T^WJ^/\vI 
For \v\ > m, with m large enough and cq large enough we have > 3 and hence 

t\v\ 



N -1> 



2y^TTW7(?iD + d/2)' 



so 



e{T) >b\[0,T]\A\ + b\A\ 

> N'^b^/TTWJ(^/\v\ + (iV + l)Db^/TT\v\^/\v\ 

> N(^Db+^bj^/TTW/(^/\v\ + ^y^l+\v\yc^ 



> (iV - 1) l^Db + -bj ^1 + I^P/cVlwl 

" 2^1 + |t;|2/c2(Z) + (i/2) V 2- 

= > /3 (> 0). 

^l + A'P/c'o{D + d/2) 

Finally we can conclude as in Case 2. 



Case 4. W^it/i a nontrivial additional force 5- 

Let us finally explain how one can take the small force ^ into account. First, we study 
the equations for \V\ and 6, where 6 is the angle between v and V{t,0,x,v). The following 
computations are valid for v large so that \V{t, 0, x, v)\ does not vanish and for a time interval 
where 6 € [-7r/2, 7r/2]. 

• For what concerns \V\, it suffices to take the scalar product with V{t, 0, x, v) of the equation 
of V. We infer 

j^\V{t,0,x,v)\^ = 2d ■V{t,0,x,v), 
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SO that 
(6.16) 



d ifj., „ x| d-V{t,0,x,v) 



Therefore, it is clear that for T > 0, there exists k > small enough, such that if the L°° 
norm of ^ is smaller than k, then one has for all {x,v) G x with \v\ > m, for any 

te[o,T], 



(6.17) 



< \V{t,0,x,v)\ < 2\v\. 



• For what concerns 6, taking the scalar product of the equation of V with v we deduce 
d 



\V{t,0,x,v)\ |7;|cos6'(t) - \V{t,0, x,v)\\v\e' (t) sme{t) 



, N N ^ (t, 0, X, • f _ 



. \V{t,0,x,v)\^ 

i H ^ 



Hence 

\V{t,0,x,v)\\v\9'{t)sme{t) 
= b(Xit, 0, x, v)) \V{t, 0, X, v) \\v\ sm{e{t)) 



V{t,0,x,v)\v\ , 
-5- ( f — —,, rr cos 0{t) . 



\V{t,0,x,v)\ 



V{t,0,x,v)\v\ 



We notice that 

- iF(t-:^:^^°^'^'^=pi^(*'°-)>^*_ 

where P{\7(t q x t;)}^ (^) denotes the orthogonal projection of v on {V{t,0,x,v)}-^ . So 



(6.18) 
Note that 



_ biX{t,0,x,v)) ^ ^ 



1 



^ \Vit,0,x,v)\^ \V{t,0,x,v)\ 



P{V{t,0,x,v)}->-('") 



\v\ sin 9{t) 



\P{V{t,o,x,v)}^i'")\ = \v\\sm{G{t))l 
so that, using Cauchy-Schwarz inequality: 

P{V{t,0,x,v)}^('") 



\V{t,0,x,v)\ 



5- 



\v\ sin 9{t) 



> 



1 



\Vit,0,x,v)\ 



mu 



Now it is rather straightforward to revisit Case 3. to include the small force ^. A helpful 
ingredient is given by the following Gronwall estimate, in which we compare the characteristics 
{X, V) associated with ^+b{x)v'^ with the characteristics {X, V) associated with the magnetic 
field b{x)v'^ alone. We have: 

\V{a + t,a,x,v) -V{a + t,a,x,v)\ < ||5||oo exp(||b||H/i.oo(l + 2|t;|)t), 
\X{a + t,a,x,v) - X{a + t,a,x,v)\ < t\\^\\ oo 6Xp(||b||^yl,cx) (1 + 2|t>|)t), 

which allows to get the result by using a perturbation argument (we recall that ||5^||oo ^ ^)- D 



(6.19) 



ON THE CONTROLLABILITY OF THE RELATIVISTIC VLASOV-MAXWELL SYSTEM 



43 



References 

1. K. Asano, On local solutions of the initial value problem for the Vlasov- Maxwell equation, Comm. Math. 
Phys. 106 (1986), no. 4, 551-568. MR 860309 (87k:76071) 

2. K. Asano and S. Ukai, On the Vlasov-Poisson limit of the Vlasov-Maxwell equation, Patterns and waves, 
Stud. Math. AppL, vol. 18, North-Holland, Amsterdam, 1986, pp. 369-383. MR 882384 (88c:35138) 

3. C. Bardos, G. Lebeau, and J. Ranch, Un exemple d'utilisation des notions de propagation pour le contrdle 
et la stabilisation de problemes hyperboliques. Rend. Sem. Mat. Univ. Politec. Torino (1988), no. Special 
Issue, 11-31 (1989), Nonlinear hyperbohc equations in applied sciences. MR 1007364 (90h:35030) 

4. , Sharp sufficient conditions for the observation, control, and stabilization of waves from the bound- 
ary, SIAM J. Control Optim. 30 (1992), no. 5, 1024-1065. MR 1178650 (94b:93067) 

5. F. Bouchut, F. Golse, and C. Pallard, Classical solutions and the Glassey- Strauss theorem for the 3D 
Vlasov-Maxwell system, Arch. Ration. Mech. Anal. 170 (2003), no. 1, 1-15. MR 2012645 (2004i:82062) 

6. N. Burq and P. Gerard, Condition necessaire et suffisante pour la controlabilite exacte des ondes, C. R. 
Acad. Sci. Paris Ser. I Math. 325 (1997), no. 7, 749-752. MR 1483711 (98j:93052) 

7. J.-M. Coron, On the controllability of the 2-D incompressible Navier-Stokes equations with the Navier slip 
boundary conditions, ESAIM Controle Optim. Gale. Var. 1 (1995/96), 35-75 (electronic). MR 1393067 
(97e:93005) 

8. , Control and nonlinearity. Mathematical Surveys and Monographs, vol. 136, American Mathemat- 
ical Society, Providence, RI, 2007. MR 2302744 (2008d:93001) 

9. P. Degond, Local existence of solutions of the Vlasov-Maxwell equations and convergence to the Vlasov- 
Poisson equations for infinite light velocity. Math. Methods Appl. Sci. 8 (1986), no. 4, 533-558. MR 870991 
(88a:76054) 

10. B. Dehman and G. Lebeau, Analysis of the HUM control operator and exact controllability for semilinear 
waves in uniform time, SIAM J. Control Optim. 48 (2009), no. 2, 521-550. 

11. R. J. DiPerna and P.-L. Lions, Global weak solutions of Vlasov-Maxwell systems, Comm. Pure Appl. Math. 
42 (1989), no. 6, 729-757. MR 1003433 (901:35236) 

12. S. Ervedoza and E. Zuazua, A systematic method for building smooth controls for smooth data. Discrete 
Gontin. Dyn. Syst. Ser. B 14 (2010), no. 4, 1375-1401. MR 2679646 

13. O. Glass, On the controllability of the Vlasov-Poisson system, J. Differential Equations 195 (2003), no. 2, 
332-379. MR 2016816 (2005a:93016) 

14. O. Glass and D. Han-Kwan, On the controllability of the Vlasov-Poisson system in the presence of external 
force fields, J. Differential Equations 252 (2012), no. 10, 5453-5491. 

15. R. Glassey and J. Schaeffer, The "two and one-half- dimensional" relativistic Vlasov Maxwell system. Comm. 
Math. Phys. 185 (1997), no. 2, 257-284. MR 1463042 (981:35143) 

16. , The relativistic Vlasov-Maxwell system in two space dimensions. I, II, Arch. Rational Mech. Anal. 

141 (1998), no. 4, 331-354, 355-374. MR 1620506 (99d:82071) 

17. R. Glassey and W. Strauss, Singularity formation in a collisionless plasma could occur only at high veloc- 
ities, Arch. Rational Mech. Anal. 92 (1986), no. 1, 59-90. MR 816621 (87j:82064) 

18. , High velocity particles in a collisionless plasma, Math. Methods Appl. Sci. 9 (1987), no. 1, 46-52. 

MR 881551 (88e:35158) 

19. R. Glassey and W. A. Strauss, Absence of shocks in an initially dilute collisionless plasma, Comm. Math. 
Phys. 113 (1987), no. 2, 191-208. MR 919231 (88k:76034) 

20. R. Glassey and Walter A. Strauss, Large velocities in the relativistic Vlasov-Maxwell equations, J. Fac. Sci. 
Univ. Tokyo Sect. lA Math. 36 (1989), no. 3, 615-627. MR 1039487 (91b:82059) 

21. S. Klainerman and G. StafRlani, A new approach to study the Vlasov-Maxwell system, Commun. Pure 
Appl. Anal. 1 (2002), no. 1, 103-125. MR 1877669 (2003a:82065) 

22. J.-L. Lions, Contrdlabilite exacte, perturbations et stabilisation de systemes distribues. Tome 1, Recherches 
en Mathematiques Appliquees [Research in Applied Mathematics], vol. 8, Masson, Paris, 1988, Controla- 
bilite exacte. [Exact controllability]. With appendices by E. Zuazua, C. Bardos, G. Lebeau and J. Ranch. 

23. G. Pallard, On the boundedness of the momentum support of solutions to the relativistic Vlasov-Maxwell 
system, Indiana Univ. Math. J. 54 (2005), no. 5, 1395-1409. MR 2177106 (20071:35231) 

24. K.-D. Phung, Contrdle et stabilisation d'ondes electromagnetiques, ESAIM Control Optim. Gale. Var. 5 
(2000), 87-137 (electronic). MR 1744608 (2001a:93051) 



44 



OLIVIER GLASS AND DANIEL HAN-KWAN 



25. J. Rauch and M. Taylor, Exponential decay of solutions to hyperbolic equations m bounded domains, Indiana 
Univ. Math. J. 24 (1974), 79-86. MR 0361461 (50 #13906) 

26. D. L. Russell, Controllability and stabilizability theory for linear partial differential equations: recent 
progress and open questions, SIAM Rev. 20 (1978), no. 4, 639-739. MR 508380 (80c:93032) 

27. J. Schaeffer, The classical limit of the relativistic Vlasov-Maxwell system, Comm. Math. Phys. 104 (1986), 
no. 3, 403-421. MR 840744 (87j:82065) 

28. , A small data theorem for collisionless plasma that includes high velocity particles, Indiana Univ. 

Math. J. 53 (2004), no. 1, 1-34. MR 2048181 (2005f:35300) 

29. S. WoUman, An existence and uniqueness theorem for the Vlasov-Maxwell system. Comm. Pure Appl. 
Math. 37 (1984), no. 4, 457-462. MR 745326 (85k:82030) 

30. , Local existence and uniqueness theory of the Vlasov-Maxwell system, J. Math. Anal. Appl. 127 

(1987), no. 1, 103-121. MR 904213 (881:35167) 

O. Glass, CEREMADE, UMR CNRS 7534, Universite Paris-Dauphine, Place du Marechal de 
Lattre de Tassigny 75775 Paris Cedex 16, France. 
E-mail address: glass@ceremade.dauphine.fr 

D. Han-Kwan, DMA, UMR CNRS 8553, Ecole Normale Superieure, 45 rue d'Ulm, 75230 Paris 
Cedex 05, France. 

E-mail address: daniel.han-kwan@ens.fr 



